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PREFACE 

During  the  past  two  years  the  author  has  been  asso- 
ciated with  Contract  N7onr-3!>8lO  between  the  Office  of  Naval 
Research  and  Brown  University,  first  as  Visiting  Professor  of 
Applied  Mathematics  at  Brown  University,  and  later  as  consul- 
tant. During  that  time,  problems  of  uniform  slabs  with  cut- 
outs and  reinforced  cutouts  have  been  studied.  The  following 
report  is  a survey  of  those  results  which  have  been  obtained 
by  the  theory  of  limit  analysis.  Some  additional  results  in- 
volving complete  solutions  [4,  8]°*^  are  not  included.  How- 
ever, all  other  results  obtained  are  either  included  in  detail 
here,  or  specific  reference  is  made  to  a previous  report. 

With  these  exceptions,  then,  the  present  survey  supersedes  the 
nB.ll''  reports  [1-9]  listed  in  the  references. 

This  survey  report  is  being  issued  at  this  time  pri- 
marily because  the  contract  supporting  it  is  about  to  expire, 
rather  than  because  the  subject  may  be  regarded  as  fully  under- 
stood. For  this  reason,  many  of  the  results  presented  are 
admittedly  incomplete.  However,  it  has  seemed  worthwhile  to 
indicate  possible  lines  for  future  attack,  rather  than  to  limit 
the  report  to  entirely  solved  problems. 

The  author  would  like  to  take  this  opportunity  to 
express  his  appreciation  to  E,  Levin  and  R.  K.  Froyd,  graduate 
students  at  the  University  of  California,  Los  Angeles,  for 


0.3  Numbers  in  parentheses  refer  to  the  references  collected 
at  the  end  of  the  report. 
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their  assistance  throughout  the  past  year  and  a half.  In  addi- 
tion to  results  for  which  they  are  directly  responsible  [5,  6, 

9 ] , they  have  helped  considerably  with  many  of  the  concepts 
and  calculations  throughout  the  report.  The  author  also  wishes 
to  thank  the  staff  at  the  Graduate  Division  of  Applied  Mathe- 
matics, in  particular  Professors  W.  Prager,  H.  J,  Weiss,  and 
P.  S,  Symonds  for  their  many  useful  suggestions  and  helpful 
criticisms  throughout  the  past  two  years. 

Philip  G,  Hodge,  Jr, 

Los  Angeles 
May,  1953 
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INTRODUCTION 

The  problems  to  be  discussed  in  this  report  are  best 
exemplified  by  an  example.  Let  us  consider  first  a thin  plane 
slab  of  thickness  h,  made  of  a homogeneous  material,  and  loaded 
in  its  plane  by  uniform  tensions  of  magnitude  2kXhTx  , 2kXhTy 
per  unit  length  applied  along  the  boundary  of  the  slab  (Fig. 
0.1a).  As  X is  slowly  increased  from  zero,  there  will  be  some 
critical  value  X^**  for  which  the  slab  will  be  no  longer  serv- 
iceable. 

We  shall  return  to  a more  precise  formulation  of  X-^** 
shortly.  First,  however,  let  us  consider  a second  slab  iden- 
tical with  the  first  except  that  it  contains  a centered  cutout 
(Fig.  0.1b).  If  this  slab  is  also  loaded  with  tensions  2kXhTx, 
2kXhTy  and  X is  slowly  increased  from  zero,  the  critical  load 
X2**  will,  in  general,  be  less  than  X^**  , since  material  which 
formerly  carried  part  of  the  load  has  been  cut  out.  Finally, 
if  thin  reinforcing  rings  are  welded  to  either  side  of  the  slab 
(Fig.  0.1c)  the  critical  load  X^*+  will  be  increased  over  the 
unreinforced  slab. 

The  design  problem  implicit  in  the  preceding  paragraph 
is  two-fold.  First,  to  design  the  reinforcing  rings  so  that  the 
slab  with  the  reinforced  cutout  has  the  same  critical  load  as 
the  original  whole  slab  (obviously  one  cannot,  in  general, 
reinforce  the  slab  to  more  than  original  full  strength).  Second, 
if  design  to  full  strength  is  impracticable,  to  determine  the 
actual  strength  of  the  final  reinforced  cutout  slab. 
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In  order  to  approach  those  problems,  it  is  first 
necessary  to  obtain  a more  precise  idea  of  the  critical  load 
X**.  To  this  end,  let  us  consider  more  closely  the  bohavior 
of  the  unrcinforced  slab  in  Fig.  0,1b.  As  X is  increased  from 
zero,  the  material  will  first  be  everywhere  elastic.  The  rela- 
tion between  the  applied  load  and  any  particular  displacement 
will  bo  linear,  and  if  the  load  is  removed  the  slab  will  return 
to  its  original  undeformed  state.  This  will  continue  until  X 
reaches  a value  X*  , at  which  point  the  slab  will  become  plastic 
at  the  most  highly  stressed  portions  A (Fig.  0.2a).  Further 


(Contained  plastic 
deformation) 


Fig.  0.2a 

Partially  plastic  slabs. 

increases  of  the  load  will  no  longer  be  rolated  linearly  to 
displacements,  and  an  unloading  would  produce  residual  stresses 
and  strains. 

For  certain  applications,  it  may  be  that  this  non- 
linear bohavior  with  residual  stresses  and  strains,  would  render 
the  slab  no  longer  useful.  However,  the  deformations  at  this 
point  are  still  small,  and  thore  is  no  tendency  for  them  to  in- 
crease rapidly.  Therefore,  it  may  be  that  further  increases  in 
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load  are  acceptable.  In  this  case,  regions  of  plastic  material 
will  begin  to  form  about  the  points  A (shaded  regions  in  Fig. 
0.2a,  0.2b),  and  possibly  to  start  at  other  points.  However, 


(Contained  plastio 
deformation) 


Fig.  0.2b 

Partially  plastio  slabs. 

there  will  still  be  a constraining  framework  of  elastic  material, 
so  that  the  slab  can  stand  further  increases  in  load.  As  X 
is  further  increased,  those  plastic  regions  will  grow,  until 
for  some  critical  value  X**  , the  elastic  framework  is  no 
longer  able  to  constrain  the  plastic  flow  (Fig.  0.2c);  even  if 


Fig.  0.2o 

Partially  plastio  slabs. 
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the  material  in  the  elastic  rogions  were  replaced  by  rigid 
material  the  overall  deformations  of  the  slab  could  still  in- 
crease to  large  values.  If  the  material  is  perfectly  plastic, 
the  plastic  regions  can  support  no  further  load,  so  that  the 
slab  as  a whole  may  deform  indefinitely.  For  a strain-harden- 
ing material  somewhat  greater  loads  can  be  supported,  but  a 
relatively  large  increase  in  deformation  will  accompany  a com- 
paratively small  increase  in  load.  Typical  load  displacement 
curves  for  both  types  of  material  are  shown  in  Fig.  0.3. 

Throughout  this  report,  we  shall  assume  that  the 

0 li 

material  is  perfectly  plastic  and  incompressible.  * For  such 
a material  we  can  specify  the  requirements  of  a ,,complete,, 
solution  of  the  problem  corresponding  to  any  value  of  the  slow- 
ly increasing  load.  Assuming  the  slab  to  be  in  a state  of 
plane  stress,  such  a solution  will  consist  of  determining  three 

stress  components  <j  , a , * and  three  velocity  components 

x y xy 

vx»  vy»  vz  a-M  as  functions  of  x and  y.  The  stress  components 
must  everywhere  satisfy  the  yield  inequality 

F(V  V V)i0  (0,1) 

where  F is  a given  function  for  the  material,  and  the  equil- 
ibrium equations 

0dx  °V  „ «TXv  ®<V  « 

+Ty  *°>Tx  + Tljr  * °*  (0.2) 


0.4  This  last  restriction  is  not  vital  to  the  general  conclu- 
sions.  but  considerably  simplifies  the  mathematical  pre- 
sentation. 
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load 


(perfectly  plastic 
material) 


load 


( Strain-hardening 
material) 


Fig.  0.3 

Load  displacement  curves. 
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Tho  velocity  components  must  satisfy  the  condition  of  incom 
pres3ibility 


Ox 


fly 


* 0. 


(0.3) 


Where  tho  strict  inequality  holds  in  0.1,  the  mater- 
ial is  said  to  be  elastic  and  must  satisfy  Hooke's  law: 


E 


E2  (J 


- 1/2  3 


y * 


E 


0v* 


<5y  - 1/2  0X 


(0.4) 


Together  with  Eqs.  0.2  and  0,3  this  provides  the  necessary  six 
equations.  Where  the  equality  holds  in  0.1  the  material  is 
said  to  be  plastic,  and  must  satisfy  the  plastic  potential 
stress  strain  law 


0P/8dx 

5i0Vx/0X)  - fdx  - t dy) 


0F/0 a v 

S(0vy/0y)  - (5y  - t 3X) 

(0.5) 

9F/drxv 

* G^v^y  + 0Vy/0x)  - TXy  * 


Equations  0.1,  0.2,  0.3  and  0,5  provide  tho  necessary  six  equa- 
tions in  the  plastic  case.  In  addition,  the  stress  components 
must  be  in  equilibrium  with  tho  appliod  loads,  and  certain 
continuity  conditions  must  be  satisfied  in  the  interior  of  the 
slab.  These  latter  conditions  determine  tho  position  of  the 
boundary  between  the  elastio  and  plastic  regions. 

A cursory  examination  of  Eqs.  0.1  through  0.5  in- 
dicates that  the  problem  of  finding  a complete  solution  is 
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extremely  difficult,  since  not  only  must  we  solve  six  equations 
in  as  many  unknowns,  but  some  of  the  equations  are  non-linear, 
and  a free  boundary  must  be  determined.  The  only  known  solu- 
tions are  in  cases  where  symmetry  or  other  restrictions  reduce 
the  number  of  independent  and  dependent  variables.  An  example 
of  such  a solution  is  given  in  (4).  In  general,  however,  Eqs. 
0.1  through  0.5  have  defied  solution,  to  date. 

Although  the  determination  of  a complete  solution 
appears  impossible,  we  can  use  recent  theorems  by  Drucker, 
Greenberg,  and  Prager  [ 10,  11  1 to  determine  something  about  the 
critical  load  X**  , A precise  formulation  of  these  theorems 
will  be  given  in  the  first  chapter.  However,  the  basic  ideas 
behind  them  may  be  stated  as  follows.  In  addition  to  consider- 
ing X**  as  the  maximum  load  which  the  slab  can  withstand  with- 
out yielding  indefinitely,  we  may  also  consider  it  as  the 
minimum  load  for  which  increased  deformation  becomes  possible 
with  no  increase  in  load.  Using  this  viewpoint,  it  can  then 
be  shown  that  by  considering  essentially  only  those  equations 
which  contain  stresses  alone  we  can  obtain  a lower  bound  on  the 
collapse  load,  while  considering  only  the  equations  which  con- 
tain velocities  we  can  obtain  an  upper  bound.  If  the  two 
bounds  are  equal,  we  have  the  actual  value  of  the  collapse  load, 
but  even  if  this  is  not  the  case  we  may  be  able  to  obtain 
bounds  sufficiently  accurate  for  engineering  purposes. 

The  body  of  this  report  will  be  concerned  with  the 
application  of  the  theorems  of  Drucker,  Greenberg,  and  Prager 
to  the  problem  of  analysis  and  design  of  slabs  with  both 
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unreinforced  and  reinforced  cutouts.  In  Chapter  I,  we  shall 
state  the  theorems  precisely  and  indicate  the  precautions 
which  must  be  used  in  applying  the  results  to  actual  situations. 
Also  we  shall  show  that  for  a symmetric  cutout  it  is  necessary 
to  examine  only  uniaxial  and  uniform  biaxial  loadings  to  de- 
duce certain  results  concerning  any  type  of  uniform  loading. 
Finally,  we  shall  discuss  certain  relationships  between  the 
reinforced  and  unreinforced  cutouts. 

In  Chapter  II  we  shall  be  concerned  with  problems  of 
narrow  reinforcing  rings  which  may  be  reasonably  approximated 
by  curved  beams.  Various  shaped  cutouts  and  types  of  rein- 
forcements will  be  considered.  Chapters  III  and  IV  are  con- 
cerned with  different  methods  of  finding  upper  and  lower  bounds, 
respectively,  for  slabs  with  unreinforced  cutouts,  while  Chap- 
ter V gives  examples  of  these  techniques  as  applied  to  various 
shaped  slabs  and  cutouts.  Finally  in  Chapter  VI  we  consider 
slabs  with  reinforced  cutouts  and  compare  some  of  the  results 
with  experimental  data. 
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I.  BASIC  CONCEPTS 

1.  Mathematical  Model.  In  connection  with  some  of 
the  results  which  will  be  presented  in  later  chapters,  the 
question  may  well  be  raised  as  to  whether  the  assumptions  made 
are  "realistic. 11  The  subject  of  what  constitutes  a realistic 
assumption  is  a most  important  one,  yet  one  which  is  often 
neglected  in  the  arguments  generated  pro  and  con  when  a new 
theory  is  presented. 

The  first  point  which  must  be  emphasized  is  that  no 
mathematical  theory  can  be  entirely  realistic.  Consider  for 
example  the  uncut  slab  in  Fig.  0.1a,  and  compare  an  actual  slab 
being  stressed  in  a testing  machine  (in  uniform  tension,  say). 
This  actual  slab  is  not  a perfect  rectangle,  since  there  will 
inevitably  be  minor  cuts  and  irregularities  along  the  edges; 
by  the  same  token  the  faces  are  not  perfect  planes.  The  applied 
load  can  be  only  approximately  uniform.  Any  mathematical 
analysis  must  moreover  assume  knowledge  concerning  the  physical 
properties  of  importance,  in  particular  in  this  case  concerning 
the  yield  stress  function  and  the  post-yield  behavior  of  the 
metal.  It  !•  veil  known,  however,  both  that  these  properties 
are  highly  complex  and  that  they  may  vary  appreciably  for 
specimens  cut  from  different  points  of  a structure,  because  of 
unavoidable  differences  in  treatment  during  the  manufacturing 
process. 

Since  it  is  practically  impossible  and  in  fact  un- 
desirable because  of  the  tremendously  complicated  nature  of 
the  analysis,  we  do  not  aven  attempt  to  give  a completely 
realistic  mathematical  account.  Instead,  we  begin  our 
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mathematical  analysis  by  postulating  a mathematical  model. 

For  the  example  considered  in  Fig.  0.1a,  the  mathematical 
model  consists  of  a perfectly  rectangular  slab  with  perfectly 
plane  faces  acted  upon  only  by  perfectly  uniform  forces  ap- 
plied precisely  at  the  edges.  In  general,  a mathematical  model 
singles  out  those  features  of  the  physical  problem  which  are 
considered  important  by  the  investigator,  replaces  them  by 
corresponding  mathematical  idealizations,  and  ignores  all  fea- 
tures considered  unessential.  Once  the  mathematical  model  has 
been  chosen,  one  operates  on  it  according  to  the  laws  of  math- 
ematics, in  order  to  predict  the  results  of  certain  experiments. 
These  experiments  are  also  carried  out  for  the  original  physical 
problem,  and  the  two  results  are  compared.  If  this  comparison 
is  satisfactory  for  the  purpose  for  which  the  results  are  in- 
tended, and  of  course  for  a sufficient  variety  of  tests,  the 
mathematical  model  is  deemed  satisfactory. 

Suppose,  however,  that  the  comparison  of  results  is 
unsatisfactory.  Even  in  this  case,  the  model  may  be  used  to 
advantage  in  any  of  several  methods.  It  may  be  that  although 
the  results  are  not  satisfactory  in  all  generality,  t-ey  are 
useful  in  a more  limited  context.  Thus,  a mathematical  model 
of  a metal  which  specifies  a linear  relation  between  stress 
and  strain  is  inadequate  for  very  large  stresses,  but  highly 
satisfactory  for  sufficiently  small  ones. 

Alternatively,  the  experience  gained  in  formulating 
and  comparing  the  model  may  suggest  ways  of  specifying  a 
better  model  and  may  suggest  the  mathematical  techniques 
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necessary  to  work  vith  the  more  complex  model.  And  finally, 
the  crude  model  may,  when  supplemented  by  suitably  chosen 
experimental  results,  enable  one  to  make  better  predications 
than  the  experiments  alone  could  lead  to. 

With  the  preceding  discussion  1*1  as  background,  we 
shall  formulate  appropriate  mathematical  models  for  slabs  with 
cutouts  in  the  following  chapters.  However,  there  are  certain 
general  features  which  are  common  to  all  cases.  First,  that 
part  of  the  slab  which  is  not  covered  by  the  reinforcement 
(and  which  we  shall  refer  to  as  the  base  slab)  is  assumed  to 
be  in  a state  of  plane  stress.  The  reinforcement  together 
with  that  part  of  the  slab  covered  by  the  reinforcement  will 
be  called  the  hub.  The  state  of  stress  within  the  hub  will  be 
differently  specified  for  different  applications,  but  in  all 
cases  the  tractions  transmitted  to  the  hub  from  the  base  slab 
will  be  assumed  to  be  uniformly  distributed  over  the  total 
thickness  of  the  hub.  Observe  that  this  assumption  may  depart 
rather  far  from  reality  if  the  reinforcements  are  thick  com- 
pared to  the  slab. 

Finally,  it  will  be  assumed  that  the  material  of  base 
slab  and  hub  is  perfectly  plastic,  incompressible,  and  satis- 
fies Tresca’s  yield  condition  [13  1 of  maximum  shearing  stress i 

F(«x,  0y,  ^xy)  » max  ft  a1  |,|<j2  I,  I *1  - c2  |]  - 2k£0  (1.1) 


1.1  For  a more  complete  treatment  of  the  subject  of  mathe- 
matical models,  see,  for  Instance,  Ref.  12,  Sec.  1.1. 
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where  and  d2  are  the  principle  stresses  in  the  slab.  A 
perfectly  plastic  material  is  one  which  behaves  elastically 
for  stresses  less  than  the  yield  stress,  flows  indefinitely 
under  constant  yield  stress,  and  cannot  support  stresses  greater 
than  the  yield  stress.  The  stress-strain  curve  for  such  a 
material  is  given  in  Fig.  1.1a  while  Fig.  1.1b  shows  a graph- 
ical representation  of  Tresca* s yield  condition.  In  Fig.  1.1b, 
states  of  elastic  stress  are  represented  by  points  within  the 
hexagon  A3CDEF,  and  the  material  cannot  support  stresses  repre- 
sented by  points  outside  of  t’  e hexagon.  Points  on  the  hexa- 
gon represent  plastic  states  of  stress*  A is  uniaxial  tension, 

B is  equal  biaxial  tension,  G is  pure  shear,  etc. 
d 


Pig.  1.1 

Perfectly  plastic  material 
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2*  The  Drucker-Oreenberg-Prager  Theorems.  As  pre- 
viously stated,  the  basic  concept  behind  the  theorems  of 
Drucker,  Greenberg,  and  Prager  is  to  determine  the  smallest 
load  for  which  an  increase  in  deformation  can  occur  with  no 
increase  in  applied  load.  T'  is  load  will  be  referred  to  as 
the  yield  load  and  the  structure  will  be  considered  safe  if 
the  acting  load  is  less  than  the  yield  load.  Subject  only  to 
the  additional  qualification  that  loads  up  to  this  value  can 
be  considered  to  be  applied  to  the  undeformed  boundary,  certain 
significant  results  have  been  obtained.  It  is  beyond  the  scope 
of  the  present  report  to  reproduce  the  proofs  of  these  theorems 5^ 
we  shall  merely  state  the  results  insofar  as  they  apply  to  t’  e 
problems  under  consideration. 

First,  then,  let  us  consider  a stress  field  (i.e.,  a 

set  of  three  functions  <L,  6 , % ) which  satisfy  the  equations 

* y xy 

of  equilibrium  (^qs.  0.1)  and  the  Tresca  yield  restriction  (1.1). 
If  any  such  stress  field  is  called  statically  admissible,  then 
the  first  theorem  may  be  stated  as  follows* 

Theorem  1.  The  yield  load  is  the  largest  load  for 
which  there  exists  a statically  admissible  stress 
field. 

It  follows  from  Thm.  1 that  any  load  for  which  a 
statically  admissible  stress  field  can  be  found  must  represent 
a lower  bound  for  the  yield  load.  Further,  it  should  be  noted 
that  no  explicit  requirements  of  continuity  were  made.  Of 


2.1  Both  theorems  were  originally  proved  for  a Prandtl-Reuss 
material  in  plane  strain  [ 10  1 • Complete  proofs  for  these 
oases  are  also  given  in  f l1*  ] , Secs.  33  and  39.  They  have 
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V> 


course,  equilibrium  implies  that  the  tractions  across  any 
element  must  be  continuous,  but  it  will  frequently  be  useful 
to  consider  stress  fields  across  which  the  remaining  stress 
component  (the  interior  component  of  stress  fl4,  p.  155  1 ) 
exhibits  a finite  jump. 

The  second  limit  analysis  theorem  is  concerned  with 
velocity  fields  which  satisfy  the  equation  of  incompressibility 
(Eq.  0.3).  Let  such  a velocity  field  be  considered  as  produc- 
ing purely  plastic  strain  rates,  and  define  a stress  field  by 
the  stress-strain  law  (Eqs.  0.5  with  the  stress  rates  equal  to 
zero).  The  resulting  internal  rate  of  dissipation  of  energy 
has  been  shown  f 1 1 to  be  given  at  each  point  by 

D « 2k  max  | e | (2.1) 

where  max  I l I denotes  the  absolutely  largest  principal  com- 
ponent of  plastic  strain  rate.  The  total  internal  energy  is 
defined  by 

<£[  **  J D dv  = 2k  J*  max  I e | dv.  (2.2) 

Also,  if  v represents  the  vector  velocities  of  points  on  the 

■<4 

boundary,  and  F the  applied  force  on  the  boundary,  the  external 
rate  of  dissipation  of  energy  is  given  by 

- |>  F • H is.  (2.3) 


2.1(Cont'd).  also  been  proved  for  the  general  perfectly  plastic 
solid  ill  J.  A detailed  account  for  the  case  of  a material 
satisfying  Tresca* s yield  condition  may  be  found  in  11  ] . 
Also,  see  [li  for  a more  thorough  discussion  as  applied  to 
states  o t plane  stress. 
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With  these  preliminaries,  an  incompressible  velocity 
field  will  be  called  kinematically  admissible  if  the  internal 
rate  of  dissipation  of  energy  does  not  exceed  the  external  rates 

A*  A*  (2<4) 

The  second  theorem  then  reads  as  follows. 

Theorem  2 . The  yield  load  is  the  smallest  load  for 
which  it  is  possible  to  find  a kinematically  admis- 
sible velocity  field. 

It  follows  from  Thm,  2 that  any  load  for  which  it  is 
possible  to  find  a kinematically  admissible  velocity  field  is 
an  upper  bound  on  the  yield  load.  Further,  it  follows  from 
Thms.  1 and  2 together  that  the  yield  load  is  the  unique  load 
for  which  it  is  possible  to  find  both  a statically  admissible 
stress  field  and  a kinematically  admissible  velocity  field.  It 
will  be  observed  that  in  this  case  also  no  explicit  requirements 
of  continuity  were  stated.  It  follows  from  the  incompressibil- 
ity of  the  material  that  the  normal  component  of  velocity  must 
be  continuous  across  any  surface,  but  the  tangential  component 
may  exhibit  a finite  Jump.  The  appropriate  internal  energy 
dissipation  may  be  evaluated  by  a limiting  process.  This  pro- 
cedure will  be  illustrated  in  Chapter  III. 

3*  Arbitrary  uniform  edge  loads.  For  the  main  part 
of  this  report  we  shall  be  considering  loadings  of  the  type 
given  in  Figs.  0.1,  l.e.,  a uniform  tensile  load  of  magnitude 
2kXhTx  per  unit  length  on  the  sides  normal  to  the  x axis,  and 
a uniform  tensile  load  of  magnitude  2kXhTy  on  the  remaining  sides. 
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If  we  consider  first  Pig.  0.1a,  and  take  X **  1 , it  is  evi- 
dent that  for  the  uncut  slab  the  resulting  stress  field  will 
be  constant i 

<JX  « 2k  Tx,  dy  « 2k  Tyt  Txy  * 0.  (3.1) 

Since  the  principle  directions  are  those  of  the  coordinate 
axes,  it  follows  from  1.1  that  the  stress  field  3.1  will  be 
statically  admissible  only  if 

max  r I Tx  |,  | Ty  |,  I Tx  - Ty  | ] ^ 1 (3.2) 

Therefore,  if  we  consider  a load  space  with  coordinates  Tx,  Ty, 
the  load  point  must  lie  in  the  hexagon  ABCDEF  (Fig.  3.1).  This 
hexagon  will  be  called  the  domain  of  safe  loads. 

Let  us  next  consider  the  same  slab  with  a cutout  (Fig. 
0.1b).  Obviously  the  domain  of  safe  loads  in  this  case  will  be 
some  region  R wholly  contained  in  the  hexagon  ABCDHF.  For 
reasons  which  will  be  elaborated  on  in  the  next  section,  we  are 
not  as  interested  in  the  actual  region  R as  we  are  in  the  largest 
hexagon  A' B'C'D'E'F' similar  and  similarly  placed  to  ABCDEF  which 
is  contained  in  R.  In  other  words,  we  seek  the  largest  number 
X such  that  if  Tx,  Ty  is  a safe  load  for  the  uncut  slab,  then 
the  loads  \TX,  XTy  are  safe  for  the  slab  with  cutout.  Such  a 
value  X will  be  termed  the  output  factor. 

With  this  definition  of  a cutout  factor,  we  can  show 
that  for  a slab  which  is  symmetric  with  respect  to  the  x and  y 
axes | only  a few  different  types  of  loading  need  be  considered 
in  order  to  establish  a lower  bound  on  X . To  this  end, 
suppose  that  a statically  admissible  stress  field  can  be 
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Fig.  3.1 

Domain  of  safe  loads. 

- slab  without  cutout 
— slab  with  cutout 
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constructed  for  the  uniaxial  tension  load  (XQ,  0 ),  that  a 
statically  admissible  field  S2  can  be  constructed  for  the 
uniaxial  tension  load  ( 0,  X ),  and  that  a statically  admis- 
sible stress  field  can  be  found  for  equal  biaxial  tensions 

( X*.  X ).  It  can  then  be  shown  that  a statically  admissible 
O'  o 

stress  field  can  be  constructed  for  any  load  point  lying  in 
the  hexagon  determined  by  XQ.  If,  in  addition,  the  slab  is 
symmetric  wit.’'’  respect  to  the  line  y * x,  then  the  existence 
of  S2  can  be  deduced  from  S^,  The  proof  of  this  result  follows 
from  the  convexity  of  the  load  domain  and  the  yield  surface. 
Details  may  be  found  in[  2]  , f 7 ],  or  ( 10  J . 

For  a reinforced  slab,  the  same  results  are  generally 
valid.  In  particular,  if  the  slab  is  to  be  reinforced  to  full 
strength,  the  dimensions  of  the  reinforcement  must  be  such  that 
>•=  1. 

4.  Relations  between  reinforced  and  unreinforoed 
cutouts.  For  certain  types  of  application,  it  is  desirable  to 
design  a reinforcement  for  a cutout  which  will  restore  the  slab 
to  full  strength,  independently  of  the  base  slab  dimensions. 

The  only  method  of  guaranteeing  this,  is  to  choose  a constant 
stress  field  in  the  base  slab,  since  such  a stress  field  will 
satisfy  the  boundary  conditions  along  any  parallel  boundary. 
Thus,  the  determination  of  a statically  admissible  stress  field 
In  the  base  slab  is  trivial,  and  it  remains  only  to  consider  the 
hub. 
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If  the  hub  is  of  constant  thickness  K,  then  it  may 
itself  be  considered  as  a plane  slab.  When  loads  2khTx,  2khTy 
per  unit  length  are  applied  to  the  base  slab,  the  same  tractions 
will  be  submitted  to  the  hub,  so  that  the  stress  on  the  hub 
boundary  is  2k(h/H)Tx,  2k(h/H)Ty.  Therefore,  the  analysis  of 
the  reinforced  cutout  slab  is  reduced  to  an  analysis  of  a plane 
unreinforced  slab. 

It  is  basically  for  this  reason  that,  in  analyzing  an 
unreinforced  slab  we  are  more  interested  in  obtaining  the  largest 
possible  load  hexagon  (Fig.  3»1)  rather  than  the  actual  domain 
of  safe  loads.  For,  once  the  value  of  K has  been  determined  for 
a slab  of  given  dimensions  and  cutout,  the  thickness  H of  a hub 
of  the  same  dimensions  and  cutout  necessary  to  restore  a slab 
of  thickness  h to  full  strength,  independently  of  the  dimen- 
sions of  the  base  slab  is 

H **  h/\  . (4.1) 

The  principles  underlying  the  above  remarks  remain 
valid  even  if  the  surfaces  of  the  hub  are  not  plane.  In  this 
case,  of  course,  other  means  must  be  used  to  analyze  the  hub, 
but  if  this  can  be  done,  Eq.  4.1  is  still  valid,  provided  that 
the  reinforcement  is  so  designed  that  h/H  = const. 
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II.  NARROW  REINFORCEMENTS 

5.  general  theory.  Consider  a cutout  of  arbitrary 
shape  in  a slab  of  thickness  h,  reinforced  by  arbitrary  narrow 
reinforcing  rings  (Fig.  5.1).  It  is  desired  to  find  a lower 
bound  for  the  cutout  factor  of  such  a slab,  valid  for  any  type 
of  uniform  loading  and  independent  of  the  dimensions  of  the  base 
slab. 

In  view  of  the  results  of  Secs.  3 and  4,  the  analysis 
of  this  problem  can  be  reduced  to  the  consideration  of  the  hub 
alone,  loaded  under  uniform  uniaxial  or  equal  biaxial  tensions. 
In  the  present  section  we  make  the  stipulation  that  the  dimen- 
sions of  a hub  are  such  that  it  may  be  reasonably  approximated 
by  a closed  curved  beam  and  that  shear  in  the  beam  may  be  neg- 
lected. Of  course,  such  an  assumption  can  be  incorporated  into 
the  mathematical  model  of  the  problem  for  any  dimensions.  How- 
ever, the  results  predicted  are  likely  to  be  reasonable  only  if 
the  reinforcement  is  reasonably  narrow,  and  if  the  thickness 
ratio  of  hub  to  slab  is  not  so  great  as  to  seriously  affect  the 
carrying  capacity  of  the  hub. 

Let  the  hub  be  referred  to  polar  coordinates  P,  0 with 
origin  at  the  center  of  symmetry,  if  any  - otherwise  at  some 
Interior  point.  Let  the  inner  radius  of  the  hub  be  denoted  by 

r-a(8)  (5.1) 

and  the  radial  thickness  by  6 (0)  so  that  the  outer  radius  is 


r » a(0)  +b  (0) 


(5.2) 
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The  height  of  the  hub  at  any  point  is  then  given  by 

H * K(!\0)  (5.3) 

We  shall  restrict  our  consideration  to  hubs  whose 
inner  and  outer  boundaries  are  sufficiently  near  to  circular, 
so  that  a "cross-section"  of  tie  hub  may  reasonably  be  considered 
as  a section  9 = const.  (Fig.  5.2a).  Within  the  framework  of 
beam  theory,  the  stresses  across  any  such  section  may  be  re- 
placed by  the  stress  resultants.  Since  shear  is  to  be  neglected, 
at  each  section  9 we  have  an  axial  force  N(9)  and  a bending 
moment  M(9),  Let  these  moments  be  taken  about  an  arbitrary 
value  r = b(0). 


(Dimensions  of  cross  section) 


Fig.  5.2a 

Cross  section  of  hub* 
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The  resultants  N and  M must  be  derived  from  the  axial 
stress  d exerted  across  the  section.  Since  this  stress  cannot 
be  greater  than  the  yield  stress,  the  resultants  will  be  sub- 
ject to  certain  inequalities.  To  determine  these  inequalities, 
consider  the  case  where  the  section  is  stressed  to  capacity. 

Thus,  in  Pig.  5*2b,  suppose  the  hub  section  to  be  stressed  to 
the  yield  limit  in  tension  forP>q  and  in  compression  for  T > rj  # 
If  the  yield  stress  is  denoted  by  2k,  the  resultant  force  and 
moment  are 


N(  0) 


ddA  « 21c 


r 


b-f  n 
-H  d r 


m a 


(5.4a) 


. ^ 


i 


(Fully  plastic  cross  section) 


Pig.  5.2b 

Cross  section  of  hub. 
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M(9)  « |d(b  -r  )dA  * 2k[  f -(b-DHdr  + fa(b-r)Hdr  ].  (5.4b) 

v v a U b+t| 

Equations  5*4  define  N and  M in  terms  of  a parameter 
1 which  may  take  on  any  value  between  a-b  and  a+b-b.  Similarly, 
if  the  regions  in  tension  and  compression  are  reversed,  one 
obtains  the  equations 

b+n  a+6 

N(9)  * 2k[  ' Hdr  + f -Hdr  1,  (5,5a) 

a «l  b+n 

_ flb+n  n a+6 

M(9)  = 2k  [ (b-DHdr  + -(b-DHdr  ] . (5.5b) 

1 a U b+n 

If  these  curves  are  sketched  in  on  N,  M plane,  the  resulting 
closed  curve  is  called  an  interaction  curve.  Stress  resultants 
corresponding  to  safe  stresses  must  lie  within  this  closed  curve. 

The  closed  hub  represents  a structure  with  one  degree 
of  indeterminancy,  so  that  the  forces  and  moments  at  any  point 
may  be  expressed  in  terms  of  the  applied  loads  and  a single 
redundant.  Let  this  latter  be  the  moment  at  section  0 « 0.  For 
simplicity  of  exposition,  we  shall  assume  that  the  hub  is  sym- 
metric with  respect  to  the  x and  y axes,  although  the  results 
are  easily  extended.  Equilibrium  of  the  first  quadrant  in  the 
y direction  then  demands  that 

N(0)  « 2kXhTy(a0  + 60)  (5.6) 

where  Sq  and  5q  are  the  values  at  0 ■ o of  a and  6 , respectively. 
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Next,  let  us  consider  equilibrium  of  an  arbitrary  seg- 
ment of  the  hub,  OABCD  (Fig.  5.3).  For  convenience  we  may  re- 
place the  uniform  loads  on  BC  and  CD  by  equipollent  concentrated 
forces 


Fx  = 2kXhTx  (a  +6  ) sin  0, 

Fy  * 2k\hTy  f (aQ  + 60)  - (a  +6)  cos  9 ] (5.7) 

acting  at  the  midpoints  of  BC  and  CD,  respectively.  Equilibrium 
in  the  direction  of  N(9)  then  yields 

N(9)  = Fx  sin  0 - Fy  cos  0 + N(0)  cos  0 

* 2kXh(a  + 6)  [ Tx  sin2  9 + Ty  cos2  9 ],  (5.Ca) 

while  moment  equilibrium  about  p * b furnishes 

M(9)  « M(0)  - (bQ  - b cos  0 ) N(0) 

♦ Fx  (b  - sin  8 + Fy  [ ?P.  * *°  +g(a*t,)c°8  9 -b  cos  8 ] 

- M(0)  + 2kXh(a46  )(b  - 3+&KTx  sin2  8 + Ty  cos2  0 ) 

- 2kXhTy  (a0  + &0)(b0  - 1°-^-?-°  ). 

(5.8b) 

The  analysis  problem  now  consists  of  the  following. 
Eliminate  tj  between  Eqs.  5A  and  between  Eqs.  5.5  to  find  the 
interaction  curve,  and  express  the  fact  that  the  point  with  co- 
ordinates N,  M must  lie  within  this  curve  by  means  of  inequalities. 


Bll-22 


28 


Substitute  Eqs,  5.8  into  these  inequalities  to  obtain  inequal- 
ities involving  M(O),0  ,X  and  known  quantities.  M(0)  is  now 
chosen  so  that  X will  be  a maximum  subject  to  the  condition 
that  these  inequalities  be  satisfied  for  all  0^  0 £ £ . Several 
examples  of  how  this  is  done,  together  with  the  application  of  these 
results  to  design  will  be  given  in  the  remaining  sections  of  this  chap 
8.  Cylindrical  rins  reinforcements T 2.  7 1 . If  the 
height  H(0)  of  the  hub  is  a constant  it  is  convenient  to  choose 


b to  be  the  center  of  the  section* 


b = a + 1/26  (6.1) 

Integrating  Eqs.  5.4  and  5.5  we  then  obtain 

N = - 2kHT],  M *=  - 2kH  -rj2),  (6.2a) 

N = 2kHi),  M = + 2kH  (.j£  - n2),  (6.2b) 


respectively.  After  eliminating  rj,  Eqs.  6.2  can  be  combined  to 
give  the  interaction  curve  in  the  form  (Fig.  6.1) 

8 k H | M j + N2  * (2kH6)2,  (6.3) 


Therefore,  safe  stress  resultants  must  satisfy 

8 k H | M J + N2  < (2kH6)2  (6.4) 

As  a first  example,  let  us  consider  a circular  cutout 

with  a circular  reinforcement,  so  that  a and  6 are  each  constant 

[ 2 ] , In  the  case  of  biaxial  tension,  T « T « 1 Eqs.  5.8 

a y 

reduce  to 
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M 


Pig.  6.1 

Domain  of  safe  stress  resultants  for  cylindrical  hub. 

N(©)  * 2kXh(a0  + 50)»  ^(©)  = M(0).  (6.5) 

Obviously  the  best  choice  of  M(0)  is  zero,  so  that  inequality  6.4 
becomes 

2kXh  (a0  + 60)  < 2kII60.  (6.6) 

Por  given  dimensions  aQ,  60,  H the  maximum  cutout  factor  for 
equal  biaxial  loading  is  thus 

X > <6’7) 

with,  of  course,  the  stipulation  that  X < 1, 

Por  uniaxial  loading  the  result  is  not  quite  so  trivial. 
In  this  case,  Tx  ■ 1,  Ty  ■ 0 so  that  Eqs.  5.8  reduce  to 


N(0)  * 2kXh  (a0  4 &0>  sin2  © 

M(©)  * M(0)  4 kXh  a0  (a0  4 60)  sin2  0 

Substituting  these  values  into  6.4  we  obtain 

8kH  j M(0)  + k\h  a0  (a0  4 60)  sin2  0 | 

+ 4 k2  X2  h2(a0  4 60)2  sin*  © £ (2kH60)2.  (6.8) 

The  left-hand  side  of  6,8,  considered  as  a function  of  sin2  ©, 
obviously  achieves  its  maximum  value  at  one  or  both  of  the  end- 
points sin2  0 = 0 or  sin2  © * 1.  Therefore,  6.8  will  be  satis- 
fied for  all  0 < © ^ J if  and  only  if 

8kH  [ M(0)  + k\haQ  (aQ  4 50)  ] + 4 k2  X2  h2  (a0  4 60>2  £ (2kH50> 
8kH  [-M(O)  - k\ha0(a0  + 60)  ] + 4 k2  X2  h2  (aD  + bQ)2  £ (2kH50)2, 
8kH  M(0)  £ (2kH60)2 

-8kH  M(0)<  (2kH60)2  (6.9) 

Inequalities  6,9  are  conveniently  treated  by  a semi- 
graphical  approach  in  [ 2 ] • However.,  in  view  of  the  more  complex 
examples  to  follow,  in  this  report,  it  is  instructive  to  solve 
them  by  a formal  mathematical  procedure.  To  this  end,  we  first 
solve  6.9  for  8kH  M(0),  obtaining 

- (2kH60)2  - 8MC2  Hha0(a0  4 60)  4 4 k2  X2  h2(a0  4 60)2  £ M(0) 

£ (2 k»0)2  - 8MC2  Hha0(ao  4 60)  - 4 k2  X2  h2(ao  4 60)2, 

(6.10a) 
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- (2kH60)2  <;  8kH  M©  <(2m0)2,  (6.10b) 

Now,  since  M(0)  is  a redundant  moment,  we  are  not  only  free  to 
choose  it  as  we  wish,  but  we  don't  really  care  what  it  is.  There- 
fore, the  only  requirement  we  make  of  Inequalities  6.10  is  that 
the  system  of  inequalities  is  not  incompatible,  but  actually 
possesses  a solution.  A necessary  and  sufficient  condition  that 
this  be  the  case  is  that  the  left-hand  side  of  each  continued 
Inequality  6.10  be  less  than  the  right-band  side.  The  resulting 
four  inequalities  are  all  independent  of  M(0)  and  may  be  written 

2 X2  h2(a0  + ft©)2  £ 2 (H50)2,  (6.11a) 

- 2(h60)2  - 2\IIha0  (aD  + 60)  + \2  h2  (a0  + &©)2  <°»  (6.11b) 
-2  (h&0)2  + 2XHha©  (a0  + 6Q)  + X2  h2  (aQ  + bQ)2  < 0,  (6.11c) 

- (H6q)2  ^ <H60)2.  (6. lid) 


Inequality  6. lid  is  an  identity,  and  since  X is  positive, 
6.11b  will  be  satisfied  whenever  6,11c  is.  Therefore,  X is  the 
largest  number  which  satisfies  6.11a  and  6,11c.  Which  of  these 
inequalities  provides  the  dominating  restriction  depends  upon  the 
ratio  b o/a0.  It  is  readily  verified  that 


bn 

X - ♦ V “ ^ * 2; 


„ \l  a02  + 2 6 o5'  -«o 

H HTaTVi’oJ  ” 


if  12  £2. 
a© 


(6.12) 
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Obviously,  if  X satisfies  6.12,  it  also  satisfies  6.7. 
Therefore,  it  follows  from  the  results  of  Sec.  3 that  X as  given 
by  6.12  is  the  cutout  factor  for  a slab  with  a circular  cutout 
reinforced  by  annular  rings  of  constant  height. 

In  particular,  if  we  want  to  design  a reinforcement  to 
restore  the  slab  to  full  strength,  we  may  set  X = 1 in  Eqs.  6.12 
and  solve  for  the  necessary  thickness  Hi 


H = h (1  + if  is  >2, 

5o  a0  - ’ 


(6.13) 


H = h 


a0  + bo 


„ ^ , if  -2  < 2. 

i aQ2  + 2 - a0  0 


Fig.  6.2  solid  curve,  shows  H/h  as  a function  of  —2  . 

ao 


As  a second  example,  consider  a square  hole  reinforced 
by  a square  reinforcement.  The  hub  to  be  analyzed  is  then 
bounded  by  two  similarly  situated  squares  of  half  side  a0  and 
a0  + &0.  The  equations  of  the  inner  and  outer  boundaries  in  the 
first  quadrant  are  given  by 


a = aQ  sec  9,  6 = 60  sec  © for  0 £ © < nA,  (6,l4a) 
a = a0  esc  ©,  6 = b0  esc  © for  n/k  £ © < n/2.  (6.l4b) 


Consider  first  the  case  of  biaxial  tension.  Due  to 
symmetry  it  is  here  sufficient  to  consider  only  the  first  octant, 

0 < © < *A.  Substituting  6,l4a  into  5.8,  and  taking  Tx  « Ty  * 1, 


we  obtain 


(6.15) 


N « 2kXh  (aQ  + &0)  sec  ©, 

M = M(0)  + kXha0  (a0  + &Q)  tan2  ©. 

The  substitution  of  Eqs.  6,1 5 into  the  yield  inequality  furnishes 
an  expression  which  is  linear  in  tan2  ©.  Therefore,  the  extreme 
values  must  come  at  the  endpoints  0 and  V*+  of  the  interval  of 
consideration.  Further,  it  is  obvious  from  6,11  that  both  M and 
N are  increasing  functions  of  ©.  Therefore,  the  yield  condition 
6.4  will  be  satisfied  for  all  values  of  © if  and  only  if 

8kII  M(hA)  + [ N(*A)  ]2  < (2k\&)2, 

- 8kH  M(0)  + [ N(0)  ]2  < (2kX&)2. 


In  view  of  Eqs,  6.15,  these  inequalities  may  be  written 
- (H&)2  + X2  h2  (a0  + 60>2  < 2H  M(0)/k, 

< (H&)2  - 2kHha0  X (a0  + &0)  - 2 X2  h2  (a0  + 6Q)2. 


A necessary  and  sufficient  condition  that  there  exist  an  M(0)  such 
that  the  above  inequality  is  true  is  that  the  left-hand  member 
be  less  than  the  right-hand.  Solving  the  resulting  inequality 
for  X we  obtain  finally 


» . H -a0  + ra02  ♦ 2 

X - 3h 


ao  *T0 


(6.15) 


Under  uniaxial  tension,  we  must  consider  both  the  first 
and  second  octants.  For  0 ^©  we  find  that  the  extreme 

values  may  occur  at  either  end  of  the  interval,  and  hence  are 
lead  to 
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2 H | M(0)  I A <H2  bQ2,  (6.17a) 

2 H | M(0)/k  + \ha0(a0  + &0)  | + 2 X2  h2(a0  + &0)2  ^ H2  &Q2. 

(6.17b) 

In  the  second  octant  M is  constant  and  N is  an  increasing  linear 
function  of  sin  9 which  leads  to  the  two  inequalities 

211  [ - MoA  - Mia0(a0  + bQ)  ] < H2  &02, 

2H  [ MoA  + Xha0(a0  + 60)  ] + X2  h2  (aQ  + 6Q)2  < H2  b2. 

(6.17c) 


Solving  each  of  the  six  inequalities  6.17  for  211  M(0)A  we  see 
that  Eqs.  6.17c  are  each  automatically  satisfied  by  one  of  the 
other  inequalities.  The  remaining  four  inequalities  obtained 
by  eliminating  M(0)  will  all  be  satisfied  if 


Since  any  X satisfying  6.18  also  satisfies  6.16,  it 
follows  that  6.18  is  a lower  bound  on  the  yield  load  for  any 
type  of  loading.  In  particular,  the  height  necessary  to  restore 
full  strength  under  these  conditions  is  obtained  by  setting 
X=  1 in  Eq.  6.18  and  solving  for  H (Fig.  6.2,  dashed  curve). 
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H = <T2  h(l  + ^2),  if  > \f2 
°o  ao  ^ 

H = 2h  - ■?.°...»  ..  ,ff  — , if  i < Vf 

if  a'o2  + 4 60  - a0  a° 

7.  Bevelled  ring  reinforcement  for  a circular  hole  T61  . 
As  an  example  of  a non-cylindrical  reinforcement,  let  us  con- 
sider the  bevelled  reinforcement  shown  in  Fig,  7,1,  Since  the 
section  of  the  hub  ABCC’B'A1  does  not  possess  any  vertical 
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line  of  symmetry,  it  is  convenient  to  take  moments  about  the 
line  BB'  which  separates  the  level  from. the  sloping  portion. 
Assuming  that  the  outermost  section  yields  in  tension,  Eqs.  5Aa 
take  on  two  different  forms  depending  upon  whether  the  neutral 
axis  is  to  the  right  or  left  of  BB'-.  Thus,  if  q is  positive  we 
have 

N * 2k  [ L q2/c  - 2Hn  - H b + (H  + h)  c/2  *}, 

(7.1a) 

M = 2k  [-  2L  rj2/3c  + H n2  - Hb2/2  - Hc2/6  - hc2/3  1 , 

while  if  q is  negative,  the  resultants  are 

N a 2k  [ - 2Hq  - H b + (H  + b)  c/2  ] , 

M = 2k[  Hq2  - Hb2/2  - Hc2/6  - hc2/3  ] . P.lb) 

It  is  convenient  to  write  Eqs.  7.1  in  terms  of  dimen- 
sionless variables  defined  by 

u » NL/2kc  H2, 

v = ML2/2kc2  H3, 

A = L(Hb  - He  + Lc/2)/c  H2 

B * L2(Hb2/2  + Hc2/2  - Lc2/3)/c  H3, 

t.  e Lq/cH  (7.2) 

Equations  7.1  then  become 

u = (2  - 2(  - A,  v «=  - 2/3  (3  + . B,  ( >0,  r.3a) 
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u = - 2£  - A,  v = K2  - B,  K < 0.  (7.3b) 

Similarly,  the  two  possibilities  which  produce  compression  at 
the  outermost  surface  lead  to 

u = - K2  + 2K  + a,  v = 2/3  K3  - K2  + B,  K > o (7.3c) 

u = 2 ( + A,  v = - 52  + B,  K < 0 (7.3d) 

The  curve  defined  by  Eqs.  7.3a  has  a cusp  at  £ = 1. 
Therefore,  it  defines  v as  a double  valued  function  of  u for 
u > A - 1,  while  no  real  v is  defined  for  u < A - 1.  However, 
since 

0 < L/H  <1 

only  one  branch  of  the  curve  has  physical  meaning.  In  the  case 

of  7.3a,  it  is  readily  verified  that  the  algebraically  greater 

branch  is  significant,  so  that,  solving  7,3a  for  v we  obtain 

v = va  = - (2/3  + A + B + u)  + 2/3  (1  + A + u)3/2.  (7.4a) 

Similarly,  only  the  lower  branch  of  Eq.  7.3c  is  significant  so 
that 

vc  = (2/3  + A + 3 - u)  - 2/3  (1  + A - u)3/2.  (7.4c) 

Finally,  Sqs.  7.3b  and  d yield  single-valued  functions 


vb  a (A  + u)2A  - B, 

(7Ab) 

'3  = - (A  - u)  A + B. 

(7.4d) 
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Next,  we  determine  u and  v in  terms  of  the  applied  load 
from  Eqs.  5.8.  If  the  slab  is  under  equal  biaxial  tension, 

Tx  ss  Ty  » l and  Eqs.  5.8  furnish 

N = 2kXh(a  + 6)  = 2kXh(a  + b + c),  M = M(0)  (7.6) 

obviously  the  best  choice  for  the  redundant  moment  is  M(0)  = 0, 
so  that  each  section  is  thus  stressed  in  pure  tension.  There- 
fore, we  must  have 

N < 2k  (Area)  a 2k  (bK  + ch  + Lc/2).  (7.?) 

Comparing  7.6  and  7.7  we  obtain 


% > bH  + ch  + Lc/2 

X ^ Wa+WJ 


(7.8  ) 


for  safety  in  biaxial  loading. 

In  uniaxial  tension  Ty  = 0,  Tx  = 1 Eqs.  5.8  furnish 


u = D X (,  v = Y + E \ c,  (7.9) 


where  in  addition  to  Eq.  7.2  we  have  introduced  the  notations 

D=Lh(a+b+  c)/c  H2, 

E * L2  h (a  + b + c)(a  + b - c)/2  c2  H3, 

Y = N(0)L2  + xL2/2kc2  H3 

C = sin2  0.  (7.10) 


itles 


The  substitution  of  7.9  into  7.5  yields  the  four  inequal- 
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-(2/3  + A + B + DXC)  + 2/3  (1  + A + D\C)3/2-  Y - EXC  < 0, 

-(2/3  + A + B - DXC)  +2/3  (1  + A - DXC)3/2+  Y + EXC  < 0, 

lA(A  + DXC)2  - B - Y - EXC  < 0, 

1A(A  - DXC)2  - B - Y + EXC  < 0.  (7.11) 

Inequalities  7.11  must  be  satisfied  for  all  values  of  C in 
0 < C i 1.  The  enforcement  of  this  condition  will  lead  to  cer- 
tain inequalities  not  involving  C.  Each  of  these  latter  can 
then  be  solved  for  the  redundant  quantity  Y.  The  elimination 
of  Y then  leads  to  various  inequalities  containing  X as  the  only 
unknown,  and  the  desired  value  of  X is  the  largest  one  satisfy- 
ing these  inequalities. 

Unlike  the  examples  in  the  previous  section,  the  functions 
of  C on  the  left-hand  side  of  7.11  may  attain  their  maximum  value 
at  an  interior  point  of  the  interval  0 ^ C < 1.  Therefore,  it 
is  not  feasible  to  obtain  a general  solution  to  7.11.  For  par- 
ticular nunerlaal  value*  of  the  parameters,  however,  a solution 
can,  of  course,  be  obtained.  For  the  particular  case 

H * 24,  h = 10,  L = 14, 

a = 24,  b = 11,  c = 20,  (7.12) 

Levin  and  Hodge  [6]  obtained  the  result  X * 0.98,  For  the 
same  example,  7.8  shows  that  X * 1 in  biaxial  tension,  so  that 
the  cutout  is  restored  to  at  least  98#  of  full  strength  under 
an  arbitrary  uniform  load. 
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8 • Symmetric  convex  reinforcement  for  a circular  cutout 

UL J. 

The  techniques  illustrated  in  the  previous  sections  can  be  used, 
in  theory,  to  find  the  yield  load  corresponding  to  any  shaped 
cutout.  However,  if  the  hub  profile  is  at  all  complex,  the 
interaction  curve  will,  in  general,  be  too  complicated  to  handle 
analytically.  Rather  than  resort  to  straight  numerical  computa- 
tions, it  may  be  desirable  to  approximate  the  interaction  curve. 

In  the  first  place,  we  observe  that  if  the  approximate 
domain  of  safe  stress  resultants  lies  entirely  within  the  true 
domain  of  safe  stress  resultants,  then  any  lower  bound  obtained 
by  the  approximate  curve  will  be  a true  lower  bound  for  the 
problem.  If  this  is  the  case,  the  approximation  will  be  called 
safe.  Since,  in  the  present  chapter,  we  are  only  concerned 
with  obtaining  lower  bounds,  It  is  most  desirable  to  know 
whether  or  not  a given  approximation  has  this  property.  We 
shall  Show  that  if  the  cross-section  of  the  hub  is  symmetric 
and  convex,  that  a simple  parabolic  approximation  will  always 
be  safe. 

To  this  end,  let  us  refer  the  hub  profile  to  Cartesian 
axes  x,  y which  are  the  axes  of  symmetry,  and  let  us  take 
moments  about  the  y axis  (Fig.  8.1),  Due  to  symmetry  we  need 
only  consider  the  first  quadrant  of  the  hub.  In  this  quadrant 
let 


y = f(x) 


(8.1) 


be  the  equation  of  the  hub.  The  boundary  may  contain  vertical 


or  horizontal  segments,  but  otherwise  we  assume  that  f is  con- 
tinuous and  momotonically  decreasing. 

If  x = C is  the  division  between  the  tensile  and  com- 
pressive regions,  Eqs.  5.4  and  5.5  may  be  replaced  by 


P K W2  • " C 

N = 4k  , f(x)dx  - f(x)dx  = 8k  f(x)dx  , 

-6/2  y0 

fiC  ) 6/2  K 

M = 4k  xf(x)dx  - xf(x)dx  = 8k  xf(x)dx  • (8.2) 
”6/2  J C o 


Due  to  symmetry,  it  is  necessary  to  consider  only  the  first 
quadrant  of  the  stress-resultant  plane,  so  that  the  absolute 
value  signs  in  8,2  may  be  dispensed  with.  Let 

*\  6/2  P 6/2 

N0  » 8k  f(x)dx,  Mq  = 8k  xf(x)dx  (8.3) 

Jo  ' 


o 
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The  true  interaction  curve  in  the  v,  p,  plane  is  then  obtained 
by  eliminating  C between  Eqs,  8.4-, 

Let  us  approximate  the  curve  8.4  by  the  parabola 

« 1 - v2.  (8.5) 

It  follows  from  Eq.  6.3  that  8.4-  reduces  to  8.5  for  a hub  with 
rectangular  section.  For  any  other  section,  we  wish  to  show 
that  the  curve  8.5  always  lies  below  the  true  curve  8,4. 

To  this  end,  let  us  fix  a value  of  C and  take  v as  de- 
fined in  terms  of  C by  Eq,  8.4a.  For  the  true  curve  H-  * is 
given  by  8.4-b,  while  for  the  approximate  curve,  it  follows  from 
8,5  that 

5 pC  ft  6/2  2 

H «1*  f «1-  [ f(x)dx  / f (x)dx  ] . (8.6) 

A u o u o 

Therefore,  if  we  define 

g(0  = (8.7) 

we  wish  to  show  that  g(C)  is  equal  or  greater  than  zero  for 

6/2. 
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g(0)  = g*  (0)  = 2(6/2)  = 0,  g*  (6/2)  < 0,  (8.8) 

and  that  g' (C)  has  not  more  than  two  roots  in  the  closed  in- 
terval 0 < C < 6/2,  Since  g(C)  is  continuously  differentiable 
it  follows  that  there  is  a stationary  point  at  C « 0 and  one 
interior  stationary  point.  This  latter  must  be  a relative  max- 
imum of  g,  from  the  last  condition  8.8.  Therefore,  since  g 
is  zero  at  the  endpoints  and  has  no  relative  minimum,  it  is 
everywhere  non-negative. 

The  preceding  result  may  be  stated  in  words  by  saying 
that  if  safe  stress  resultants  can  be  found  for  a rectangular 
section  they  can  be  found  for  any  other  hub  section  with  the 
same  6,  N0  and  M0.  Therefore,  a safe  approximation  for  any 
symmetric  convex  section  can  always  be  found  in  terms  of  a 
rectangular  section  and  the  analysis  completed  as  in  Sec.  6. 

In  particular,  it  is  easily  shown  [9  3,  that  under  equal 
biaxial  tensions  2kXh  per  unit  length,  we  must  have 

^ S A/C  h(a  + 6)],  (8.9) 

where  A is  the  total  area  of  the  section,  and  h the  thickness 
of  the  base  slab.  Under  uniaxial  tension  2kXh  we  must  have 

X < min  (X^,  X2),  (8.10) 

where  X^  is  given  by  8.9  and 


8.1  Details  may  be  found  in[  93 
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*2  5 »1  t'f  2 + (aN0AM0)2  - (aN0/4M0)  ],  (8.11) 

Finally,  it  follows  from  the  analysis  of  Sec.  3 that  a lower 
bound  on  the  cutout  factor  is  given  by  8.10. 

In  [9]  the  above  analysis  is  used  to  discuss  a toroidal- 
type  reinforcement  for  a circular  cutout.  In  that  example  It 
is  shown  that  the  error  introduced  by  approximating  the  inter- 
action is  negligible.  In  general,  however,  no  such  conclusion 
has  been  obtained. 

9 • Design  of  narrow  reinforcement  for  full  plasticity . 

If  the  reinforced  cutout  is  to  be  subjected  to  only  one  type 
of  loading,  then  it  may  bo  possible  to  design  a cylindrical 
reinforcement  for  the  cutout  so  that  it  will,  like  the  "one 
horse  shay"  reach  the  yield  stress  at  all  points  simultaneously. 

In  other  words,  for  a given  a(0)  wo  wish  to  choose  the  rein- 
forcement 6(G)  so  that  the  Equality  6.3  will  hold  at  all  points 
of  the  hub,  rather  than  the  Inequality  6.4. 

Since  the  general  theory  becomes  quite  complicated  [ 7 ] , 

we  shall  here  illustrate  It  by  considering  a circular  cutout 

0 1 

stressed  in  uniaxial  tension,  Further  we  wish  to  design 
the  reinforcement  so  as  to  restore  the  slab  to  full  strength. 

Thus  Tx  8 1,  Ty  * 0,  \ = 1 in  Eqs.  5.8,  so  that  the  resultants 
are 

N(Q)  = 2kh(a0  + 6)  sin2  0, 

M(Q)  * M(0)  + kha0(a0  + 6)  sin2  6.  (9.1) 

9.1  The  basic  idea  for  this  example  was  suggested  by  J.M.  English 

c 15] . 
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We  shall  verify  a posteriori  that  tho  quantity  6 sin2  © is 
an  increasing  function  of  ©,  so  that  N and  M are  both  a max- 
imum at  Q = V2.  Therefore,  if  6,3  is  to  be  everywhere  valid, 
M(0)  must  bo  everywhere  loss  than  zero  except  at  © = n/2  where 
it  will  equal  zero,  i.e.,  the  cutout  will  yield  in  pure  tension 
there.  With  M(0)  negative,  Eq.  6.3  may  be  written 

- 2IJ  M/k  + (N/2k)2  *s  H2  62.  (9.2) 

At  © = 0 Eqs.  9.1  and  9.2  show  that 

N(0)  = 0, 

2H  M(0)/k  = - H2  602,  (9.4) 

while  at  © = V2  they  furnish 

N(  n/2  )/2k  * h (a  + 8j_)  E 

2K  M(V2)A  = - H2  6q2  + 2Hha0  (a0  + 63.)  = 0,  (9.5) 

whore  6^  is  the  value  of  6 at  © = -n/2 . Eliminating  6^,  between 
Eqs.  9.5  we  obtain 


Finally,  the  substitution  of  9.1,  9.4,  and  9.6  into 
Eq.  9.2  shows  that 

[(1  + 2a02)2  - sin*  ©]  p2  + 2a0  sin2  © [ (1  + 2aQ2)  - sin2  ©]  p 
- [(1  + 2a02)2  - 2 a^2  sin20(l  + 2aQ2)  + a02  sin4  © ] = 0,  (9.7) 
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where  we  have  used  the  dimensionless  variables 


p (0)  - 6(0)/  60,  Oq  = a0A 0. 


Introducing  the  further  abbreviation 

f(0)  = sin2  0 / (1  + 2a02), 
we  may  write  the  solution  of  Eq.  9.7  as 


P = 


% f 
1 + f 


■\ 


1-2  % f/a  + f) 


i - 


(9.8a) 


(9.8b) 


(9.9) 


where  the  positive  root  is  chosen  since  p must  be  positive. 

In  Figs.  9.1  and  9.2  we  have  sketched  curves  showing  the 
shape  of  the  reinforcement  for  the  cases  6Q  = 1 (H  = 3.00  h) 
and  60  = 2 (H  = l.|>0  h)  respectively.  In  each  figure,  the  dot- 
ted curve  shows  the  circular  reinforcement  to  full  strength  for 
the  same  value  of  H,  computed  from  Eq.  6.13.  Finally,  Fig.  9.3 
shows  similar  curves  for  a square  cutout  under  various  loadings. 
The  details  for  this  latter  case  may  be  found  in  [ 7 3 . 


h.a-kfi 


1+9 


Fig,  9.1 

Full  reinforcement  for  circular  cutout.  a 1. 

1 o 
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III.  UPPER  BOUNDS  FOR  WIDE  SLABS 

10,  Discontinuous  velocity  fields.  The  results 
presented  in  the  previous  chapter  were  concerned  with  a narrow 
reinforcement  so  that  the  hub  might  reasonably  be  approximated 
by  a curved  beam.  Actually,  the  range  of  application  of  such 
techniques  is  rather  limited,  since  if  the  ratio  H/h  becomes 
at  all  large  there  is  serious  question  as  to  the  carrying 
capacity  of  the  hub,  while  if  the  ratio  6 /a  becomes  large  the 
hub  is  not  well  approximated  by  a beam. 

A more  realistic  analysis  in  the  case  of  wide  unrein- 
forced slabs  or  slabs  with  wide  reinforcements  is  to  approx- 
imate the  actual  stress  state  by  plane  stress.  In  this  and 
the  two  succeeding  chapters,  we  shall  consider  the  problem  of 
finding  the  yield  load  of  unreinforced  slabs,  assuming  a state 
of  plane  stress.  V/e  shall  use  both  of  the  Drucker-Greenberg- 
Prager  theorems  in  order  to  obtain  both  lower  and  upper  bounds. 
The  accuracy  of  the  answer  will  be  determined  by  the  relative 
closeness  of  these  two.  In  cases  where  the  lower  bound  is  far 
below  the  upper  bound,  the  problem  must  be  regarded  as  still 
unsolved. 

A very  simple  type  of  velocity  field  which  is  quite  use- 
ful is  to  assume  that  the  slab  consists  of  two  rigid  parts, 
and  that  motion  occurs  by  the  relative  sliding  of  these  two 
parts.  In  the  present  section  we  shall  develop  a general 
formula  for  computing  the  internal  rate  of  dissipation  of 
energy  for  such  a motion  (see  also  [16]), 
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We  consider  two  sections  of  material  separated  by  a 
common  plane  surface  of  area  A,  with  one  surface  in  motion 
relative  to  the  other  with  a tangential  velocity  A vt  (Fig. 
10.1a).  As  mentioned  in  Sec.  2,  such  a discontinuous  velocity 
field  must  be  considered  as  the  limit  of  a suitably  chosen 
continuous  velocity  field.  Let  us  choose  coordinate  axes  so 
that  the  x axis  is  normal  to  the  plane  of  sliding,  and  the  y 
axis  is  in  the  direction  of  sliding,  and  consider  the  con- 
tinuous motion  defined  in  Fig.  10.1b.  As  in  Fig.  10.1a, 
region  1 is  at  rest  while  region  3 moves  upwards  with  velocity 
AVj..  In  the  transition  region  2,  we  define  the  velocity  field 

vx  = vz  = 0,  v = (x/6)  Avt.  (10.1) 

v 

Observe  that  this  field  is  incompressible,  and  -chat  the  com- 
bined velocity  field  is  continuous  in  region  1,  2,  and  3. 

Since  the  velocity  fields  in  regions  1 and  3 are  each 
constant,  the  strains  are  all  zero.  Therefore,  energy  is 
dissipated  only  in  region  2.  In  region  2,  the  absolutely 
largest  principal  strain  rate  has  the  constant  value  Avt/25, 
Therefore,  3q.  2.2  yields 

* 21c  I ( Avt/26)dv  = k ( Avt/5)  V2  = kAAvt  (10.2) 

Jv2 

Independently  of  the  value  of  6 . Therefore,  it  follows  that 
Eq.  10.2  must  also  be  valid  for  the  limiting  discontinuous 
case,  5=0. 
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11.  Sliding  out  of  plane  f 2 1.  A simple  upper  bound 
for  any  plane  slab  may  be  found  immediately  by  considering  a 
shear  out  of  the  plane  at  45°.  Consider,  for  example,  such 
a shear  across  the  x axis  of  an  unreinforced  slab  (Fig.  11.1). 
Since  the  plane  of  sliding  is  at  45°  to  the  plane  of  its 
height  will  be  h \JTT  Further,  since  the  sliding  takes  place 
along  the  weakest  section  AB,  the  sliding  plane  consists  of  two 
rectangles,  each  of  length  1 - aQ.  Therefore,  it  follows  from 
Eq.  10.2  that  the  total  internal  rate  of  dissipation  of  energy 
is 

<0±  = k 2hrz'(l  - aG)  V (11.1) 

The  external  rate  of  dissipation  of  energy  is  computed 
from  Eq.  2.3.  Since  there  is  no  motion  in  the  x direction,  the 
loads  Tx  do  no  external  work.  The  projection  of  the  velocity 
vector  in  the  direction  of  the  load  T is  v/i/2,  so  that  in 
view  of  2.3, 

* | 2 k\Ty  (v/\/T)ds  * 2 VT  k h X Ty  V.  (11.2) 

According  to  Theorem  2,  an  upper  bound  for  the  load  X is  ob- 
tained if  Oe  and  are  equal.  Thus 

XTyil-a0,  (11.3) 

Since  the  average  stress  across  AB  is  2kXTy/(l  - a0),  Eq.  11.3 
shows  that  this  stress  cannot  exceed  the  yield  stress.  This  is 
in  marked  contrast  to  the  yielding  of  notched  bars  under  con- 
ditions of  plane  strain  [ 14,  Sec.  33  3. 


Fig.  11.1 

Sliding  out  of  plane. 
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12.  Sliding  In  plane.  The  concept  of  rigid  parts 
sliding  along  planes  can  also  be  used  for  the  case  where  the 
plane  of  sliding  is  perpendicular  to  the  slab.  As  an  example, 
consider  an  annular  slab  of  inner  radius  a and  outer  radius  1, 
sliding  along  four  symmetrically  placed  lines  in  the  plane,  as 
indicated  in  Fig.  12.1,  The  slab  is  under  uniaxial  tension  in 
the  y direction*  Tx  = 0,  Ty  **  1. 

Due  to  symmetry,  we  need  consider  only  the  first  quadrant 
of  the  slab.  In  this  quadrant  the  x and  y components  of  velocity 
are 

Vx  = - A,  Vy  = 0 in  region  1, 

Vx  = 0,  Vy  as  B in  region  2.  (12.1) 

However,  it  is  more  convenient  to  refer  the  motion  to  the  dir- 
ections n and  t,  normal  and  tangential  to  the  line  of  discon- 
tinuity AB  in  Fig.  12.1.  Equations  12,1  must  then  be  replaced 
by 

Vm  = - A sin  Y,  V^  = - A cos  Y in  region  1, 

Vm  as  - B cos  Y>  vt  = B sin  Y in  region  2,  (12.2) 

The  constants  A and  B in  Eqs.  12.2  are  not  independent,  since 
the  normal  components  of  velocity  must  be  the  same  in  the  two 
regions.  Thus 


3 **  A tan  Y, 

and  the  discontinuity  in  tangential  velocity  is  readily  computed 
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to  be 


A Vt  = A sea  y 


(12.3) 


Substituting  12.3  into  10,2,  we  see  that  the  internal  rate  of 
dissipation  of  energy  is 


= k h L A sec  y , 


(12.4) 


where  L is  the  length  of  AB. 

The  external  rate  of  dissipation  of  energy  is  given  by 

nV2 

,£)_  * 2 k X sin  ©Bhd©  = 2kXhA  tan  y cos  B.  (12,5) 

e up 

An  upper  bound  on  X is  then  obtained  by  equating  12.4  and  12. 6 i 


- 2 sin  y sin  P* 


(12.6) 


Note  that  both  Eqs.  12.4  and  12.6  represent  one-quarter  of  the 
total  energy  of  the  slab  so  that  they  may  be  equated. 

The  line  of  sliding  AB  has  two  free  parameters  which  are 
conveniently  taken  to  be  the  angles  a and  P as  Indicated  in 
Fig.  12,1.  In  terras  of  these  angles  we  have 


L = ^1  - 2 a cos  ( p - a)  + aP~  , 
sin  y=  (sin  P - a sin  a)/L. 


(12.7) 


Substituting  12.7  into  12,6  we  obtain 


x ✓ 1 - 2 a cos 

- 2 cos  B (sin 


3 - a)  + a 
- a sin  a ) • 


(12.8) 
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Any  values  of  a and  9,  subject  to  certain  obvious  geome- 
tric restrictions,  will  furnish  an  upper  bound  for  X.  In 
particular,  if  AB  is  assumed  to  be  a radial  line,  a = £ and 
the  right-hand  side  of  12.8  becomes  (l-a)/sin  2 {3  , Since  this 
can  never  be  less  than  1-a,  sliding  along  a radial  line  can- 
not furnish  a better  upper  bound  then  sliding  out  of  the  plane 
(Sec.  11). 

Obviously  the  best  upper  bound  of  the  type  considered 
is  obtained  by  choosing  a and  p so  as  to  minimize  the  right- 
hand  side  of  12.8.  To  this  end,  let 


If  f (u,v) 


u = 1 - 2a  cos(  P - a ) + a^ 

v = cos  P (sin  P - a sin  a ) . 

= is  to  be  a minimum,  then 
2v 

§ = 1/2  - U 8.V/9?  = 0i 


» 


0 


so  that 


(12.9) 


u _ da/da  _ 9u/g6 
v “ dv/da  ~ V/&P  * 


(12.10) 


or 


2 a cos(  B - q ) + a2 
p cos  p - a "sin  a aos^ 


= -2-&  .sin ..(  p -a  J 

cos  p cos  a cos  2p  + a sin  a sin  p * 

(12.11) 


Equating  the  last  two  and  first  two  members  of  the  continued 
equality  12,11,  y/e  obtain 
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a cos  p cos  a - a sin  {3  sin  a = cos  2 p, 

(a2  + cos  2 P)  cos  a + sin  2 (3  sin  a==  2 a cos  P.  (12.12) 


Equations  12.11  may  be  formally  identified  with  simultaneous 
linear  algebraic  equations  in  cos  a and  sin  a and  solved  to 
yield 


cos  a = 2 cos  P (cos  2 j3  + 

a (2  cos  2 p + 1 + a2)  * 

sin  a = ?.  .--a! . 

-a  sin  P(2  cos  2 p+  1 + ad) 


(12.13) 


Substituting  12.13  into  the  identity  sin^  a + cos  a = 1 
and  introducing  the  notation 


C * cos  2 p 


(12.14) 


we  obtain,  after  some  simplification, 

(1  - a2)  (2C  2 + a2  C - a2)  = 0. 


Since  a must  be  actually  less  than  one,  the  second  factor  must 
vanish,  so  that 

C=  (aA)  (-a  ± f a2'  +”8“).  (12.15) 


Now,  it  follows  from  Eq.  10  that  the  minimum  value  of 

f is 

f * 1/2  a * 1/2  = tan  p - tan  a.  (12.16) 

v ov/oa 

Replacing  a and  p in  Eq.  12.16  by  their  values  in  terms  of  £ 
as  given  by  Eqs.  12.14  and  12.13,  we  obtain  finally 
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X < f 


C tl  - a2) 

(t  + az)sr“rx 


2 


(12.17) 


where  C Is  given  in  terms  of  a by  Eq.  12.15.  It  is  readily 
verified  that  the  upper  sign  in  12.15  will  furnish  a lower 
value  of  f and  hence  is  the  correct  choice.  In  Fig.  12.2  we 
have  sketched  the  upper  bound  on  X as  a function  of  a.  Before 
accepting  this  result,  it  was  necessary  to  find  p and  a from 
Eqs.  12, 14  and  12.13  and  verify  that  they  were  geometrically 
admissible  angles.  This  was  found  to  be  the  case  for  all 
values  of  a.  Observe  that  for  all  values  of  a we  have  obtained 
an  improved  result  for  this  case  over  that  obtained  in  Sec.  11, 
as  shown  by  the  dashed  line  in  Fig.  12.2, 

13.  Bending  in  Plane.  A totally  different  type  of 
yielding  is  suggested  by  the  approximate  analysis  of  the  pre- 
ceding chapter.  If  the  slab  were  to  be  replaced  by  a curved 
beam,  the  beam  would  fail  by  means  of  yield  hinges.  That  is, 
there  would  be  certain  sections  of  the  beam  where  there  would 
be  a concentrated  rotation  and  extension  [17]  (or,  equivalently, 
an  off  center  rotation  [l8]  ),  with  the  remainder  of  the  beam 
remaining  rigid.  There  must  be  just  enough  of  these  yield 
hinges  to  form  a mechanism. 

Let  us,  then,  consider  such  a yield  hinge  from  a plane 
stress  viewpoint.  Choosing  the  y axis  to  be  normal  to  the 
beam  at  the  yield  hinge,  and  the  plane  of  bending  to  be  the 
x,y  plane,  consider  a rotation  about  the  point  A in  Fig.  13.1. 
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X 


Fig.  12.2 

Upper  bound  on  yield  load  of  annular  slab* 

_ _____  sliding  in  plane 

— — — — sliding  out  of  plane 
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Fig.  13.1 

"Yield  hinge"  in  plane  stress. 

The  rigid  portion  of  the  beam,  region  3 will  rotate  about  the 
point  A with  an  angular  velocity  ©,  Therefore 

Uj  = ■ fi  (y  - t),  Vj  = 5 x (13.1) 

where  u^  and  are  the  velocity  components  in  the  direction 
of  the  x and  y axes,  respectively.  In  region  1 u must  equal 
zero  along  AB  by  symmetry,  and  must  be  continuous  with  u^  along 
AC.  Similarly,  v must  be  continuous  along  AC,  and  there  is 
to  be  no  shearing  strain  in  the  vertical  direction.  Such  a 
velocity  field  is  given  by 

u^  ■ - © x cot  a , Vj  » i (y  4 ) tan  a.  (13.3) 
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Similar  considerations  show  that 

u2  = 6 x tan  P , v2  = © (-t  - y)  tan  {3  . (13.3) 

There  are  no  strains  in  region  3»  while  in  region  1 the 
principal  strain  rates  are 

e = - 0 cot  a,  £ = © tan  a , e = © (cot  a - tan  a ),  (13.4) 

Observe  that  the  axial  velocity  w is  thus  not  strictly  con- 
tinuous. However,  for  a sufficiently  thin  slab  its  action  is 
reasonably  neglected.  Separate  analyses  must  be  made  for 
greater  or  loss  than  For  definiteness  we  shall  take 

a < vA-  so  that  the  absolutely  largest  principle  strain  rate 
has  the  value 


max  | e | = © cot  a . (13.5) 

Substituting  13.5  into  2.2,  we  find  that  the  internal  rate  of 
dissipation  of  energy  is 

2 k © h cot  a ky  (13.6a) 

where  is  the  area  of  region  1.  Similarly,  if  {3  <flA, 

2 k © h cot  {3  A2.  (13.6b) 

Since  = 0 the  total  internal  energy  is  the  sum  of  the 

two  Eqs.  13.6. 

As  an  example,  consider  a square  hole  in  a square  slab 
under  uniaxial  tension  (Fig.  13.2).  The  velocity  of  the  rigid 
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Fig.  13.2 

"Bending"  velocity  field  for  square  slab  with  square  cutout. 


part  of  the  material  is 

u = 0 (y  - a - c),  v=-0(x-a-  f).  (13.7) 


In  this  case,  the  area  of  region  1 is 
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= 1/2  b2  tan  a , 

so  that 

= k h i b2,  (13.8) 

independently  of  the  angle  a . 

Since  we  seek,  the  lowest  possible  upper  bound,  we 
generally  desire  to  choose  the  parameters  of  a problem  so  as 
to  minimize  the  internal  and  maximize  the  external  energy  dis- 
sipation rates.  Since  is  independent  of  a , we  may 

therefore  choose  a so  as  to  maximized.  In  view  of  Eq.  2,3, 

© 

i»  1 1 1 

/)=!  2k\hvdx  = 2k\h  vdx  (13.9) 

(!  0 to 

so  that  we  wish  to  choose  v as  large  as  possible.  Since 

v-j_  = - 0 [ tan  a (y  - a - c)  - (a  + f ) ] 

is  less  than  v^  along  the  top  surface,  it  is  advantageous  to 
make  region  1 as  small  as  possible,  i.e,,  to  choose  the  limit- 
ing case  a = 0.  Therefore,  we  may  use  13.7  in  the  entire 
integral  13.9  to  obtain 

=2kXh9(a+f-  1/2).  (13.10) 

The  total  internal  energy  dissipation  rate  is  given  by  four 
terms  similar  to  13.8 1 

= k h © (b2  + c2  + d2  + f2) 

= k h 0 [b2  + (1  - a - b)2  + f2  + (1  - a - f)2  ] . (13.11) 
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Therefore  an  upper  bound  is  given  by 


X < b2  + (1  - a - b)2  + f + q - a - f)2 
- 2 a + 2 f - 1 


(13.12) 


for  any  geometrically  admissible  choice  of  b and  f. 

Since  the  denominator  is  independent  of  b,  the  best 
choice  is 


b » I.-  ft 

2 

which  locates  the  hinge  at  the  midpoint  of  the  section.  There- 
fore, 13.12  becomes 

x < x0  - 'HV  3 ' fj!  (13.13a) 

If  a < l/3j  the  minimum  value  of  X0  occurs  for  f = 1 - a,  i.e., 
the  hinge  is  at  the  outer  edge.  However,  for  1/3  < a,  \0  has 
a relative  minimum  at 

f = 1/2(1  + \[  2 - 4a  + 3a^  - a.  (13.13b) 

With  this  value  of  f,  the  resulting  XQ  represents  the  best 
upper  bound  of  the  type  considered.  Figure  13.3  shows  that  we 
have  obtained  an  improvement  over  the  result  of  Sec.  11  for 
all  a > 1/3. 


• 2 

0 0.2  0.4 

Fig,  13.3 

Upeer  bounds  for  square  slab  with  square  cutout 

___________  bending 

— — sliding  out  of  plane 
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IV.  LOWER  BOUNDS  FOR  WIDE  SLABS 
14*  Discontinuous  stress  fields  [3l*  As  we  have  seen  In  the 
previous  chapter,  it  is  not  difficult  to  construct  upper  bounds 
for  the  yielding  of  slabs  with  cutouts*  However,  the  determina- 
tion of  statically  admissible  stress  fields  is,  in  general,  more 
difficult. 

In  the  case  where  we  are  concerned  primarily  with  uni- 
form loads  on  rectangular  boundaries,  it  is  often  possible  to 
construct  a discontinuous  stress  field  consisting  of  several 
regions  of  constant  stress  separated  by  straight  lines  of  dis- 
continuity* This  approach  has  also  been  used  for  problems  in 
plane  strain  [14,  Sec*  33]* 

We  shall  find  it  convenient  to  define  the  stress  field 
in  terms  of  the  principal  stresses  and  a representative  angle* 

To  this  end,  let  0 be  the  smallest  non-negative  angle  between 
the  x axis  and  a principal  direction,  let  s be  the  principal 
stress  across  the  element  inclined  at  0 to  the  x axis,  and  let 
r be  the  other  principle  stress  in  the  plane  of  the  slab*  Let 
w be  the  reduced  mean  normal  stress  and  x the  reduced  principal 
stress  difference: 

w = (r  + s]^k,  X 55  (*  - sjfek,  (14*1) 

Let  N and  T be  the  normal  and  tangential  stresses 
across  an  element  inclined  to  the  x axis  at  an  angle  a*  In  terms 
of  the  new  variables 

= u - x cos  2(0  - a),  ^ = - x sin  2(0  - 0),  (14. 2) 

Across  any  line  of  discontinuity,  these  exterior  components  of 
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stress  must  be  continuous,  but  the  remaining  interior  component 
of  stress  may  be  discontinuous*  Thus,  if  is  the  angle 
between  the  x axis  and  the  straight  line  of  discontinuity  sepa- 
rating regions  k and  j,  we  must  have 

“j  ' Xj  cos  2(6J  - \])  = “k  - xk  2<®k  - V’ 

(14.3) 

Xj  ain  2{#J  - V * Xk  2<»k  - «kj)* 

Any  stress  field  which  satisfies  14*3  and  the  stress 
boundary  conditions  will  be  in  equilibrium.  If  it  is  to  be 
statically  admissible,  it  must  also  satisfy  the  yield  condition. 
In  terms  of  the  present  variables  this  may  be  written 

max[|JL|,  |^|,  1x13  <!•  (14.4) 

As  an  example,  let  us  consider  a square  slab  with  a 
square  cutout  under  uniform  uniaxial  tension.  The  assumed  dis- 
continuous stress  field  is  shown  in  Fig.  14. 1.  Symmetry  demands 
that  the  principal  directions  in  regions  1,  3,  and  4 be  parallel 
to  the  coordinate  axes,  so  that 

Qi  = ©3  = = 0.  (14. 5a) 

Due  to  this  symmetry  we  need  consider  only  the  first  quadrant  of 
the  slab#  In  regions  1,  3»  and  4 the  stress  normal  to  the  bound- 
ary must  be  in  equilibrium  with  the  load,  so  that 

sx  = 2k\,  s^  = 0,  r3  = °»  (14. 5b) 

There  will  be  three  pairs  of  Eqs.  14,3  corresponding 
to  the  three  lines  of  discontinuity  in  the  first  quadrant.  Since 
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Fig,  14,1 

Discontinuous  stress  fields. 
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a = %/b  these  may  be  written 
23 

r-, 

w2  * x2  cos  2^®3  " °12^  = ^2k^1  ” cos  2ai2^  + ^(1  + 003  2al2^» 
u2  * X2  c°s  2<93  - a2i+)  = <l£>(l  - cos  20^), 

S-5 

W2  - x2  sin  2e2  = 2k»  (l'+.6) 

P 

X2  sin  2(Q2  - a12)  = (X  - ^i) sin  20^, 

X2  sin  2(e2  - a2h)  = - (^)sln  20^, 

- x2  °os  2©2  = • 


The  geometry  of  the  stress  field  is  entirely  determined 
by  the  choice  of  the  parameter  x (Fig*  14-,  1),  Indeed,  we  have 


tan  a = 1 - x,  tan  o . = 1 - 


(1 4.7) 


12  ~ 24-  a 

Therefore,  Eqs.  14-.6  may  be  solved  for  the  six  unknown  quantities 
in  terms  of  x: 


tan  290  = 

2 X 7 


\/x2  + 4a2 
a x 9 

V-fh 


Xa 

(1  - a)  (1 


Xa 

2k  ~ (1  - a) ( 1 - x)  9 


— • 

a 


(14*8) 
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As  a check  on  the  computations  it  may  be  verified  that  the  result 
ing  stresses  satisfy  overall  equilibrium  requirements  across  the 
x and  y axes. 

The  problem,  then,  is  to  determine  the  largest  value 
of  X for  which  14.4  will  be  satisfied  in  each  of  the  four  regions 
In  each  region,  that  inequality  which  would  be  violated  for  the 
smallest  value  of  X will  be  called  the  governing  inequality.  In 
the  present  example  it  may  be  determined  by  inspection  in  each 
region.  The  four  governing  inequalities  obtained  by  substituting 
14, 8 into  14,4  may  then  be  written 


Xa  < (1  - a)(l  • x), 

(14.9a) 

X Vx2  + 4a2  <.  x(l  - a), 

(14. 9b) 

X < 1 - a, 

(14. 9c) 

Xa  <:  (1  - a)(x  - a). 

(14. 9d) 

\ / 2 2 

Since  x/Vx  + 4a  is  always  less  than  1,  Inequality  14. 

9c  is 

obviously  satisfied  whenever  14.9b  is,  so  that  14,9  are 

equiva- 

lent  to 

X < (1  - a)(l  - x)/a, 

(14. 10a) 

X ^xd  - a)/\/x  + 4a  , 

(14. 10b) 

X < (1  - a)(x  - a)/a. 

(14. 10c) 

Subject  to  the  geometric  restriction  a < x ^ 1,  we  may  choose  x 
in  14*10  so  as  to  obtain  the  maximum  possible  value  of  X*  While 
x may  be  formally  eliminated  from  the  inequalities  as  in  Sec,  6, 
the  resulting  functions  of  X cannot  be  handled  in  closed  form* 
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Therefore,  the  following  procedure  seems  more  appropriate. 

The  general  character  of  the  right-hand  sides  of  14.10, 
considered  as  functions  of  x in  the  interval  a < x < 1 is  inde- 
pendent of  a.  Thus  14, 10a  is  a straight  line  falling  from 

p 

(1  - a)  /a  at  x = a to  0 at  x = 1,  while  14. 10c  is  a straight 
line  rising  from  0 to  (1  - a)  /a  in  the  same  interval.  Finally, 
14, 10b  is  monotonically  increasing  and  positive  in  a < x < 1. 

These  curves  are  sketched  in  Fig.  14, 2.  We  must  distinguish  two 
cases  depending  upon  whether  curve  (b)  passes  above  or  below  the 
intersection  of  (a)  and  (c).  In  the  former  case,  we  must  take  x 
as  the  intersection  of  (a)  and  (c)  in  order  to  maximize  X,  while 
in  the  latter  case  x is  chosen  as  the  intersection  of  (b)  and  (c). 
After  some  computation  [3]  we  find 

x = I (1  + a),  X — ( 1 — a)2/2a,  if  a > 0.443, 

(14. 11) 

, X = .U.r-ftlU.i.xI,  if  a < 0.443. 

The  resulting  lower  bound  on  the  collapse  load  is  shown  in  Fig. 

14. 3»  Also  shown  is  the  result  for  biaxial  tension  which  will  be 
obtained  in  Sec.  18# 

15.  Stress  functions.  While  the  method  of  discontinuous  constant 
stress  fields  is  easy  to  apply  in  a few  examples,  it  does  not 
adapt  Itself  to  non-rectangular  domains.  A second  technique, 
which  holds  great  theoretical  promise  although  an  efficient  means 
of  exploitation  is  yet  to  be  found,  is  to  construct  a stress 
function.  As  is  well  known,  if  the  stress  components  are  derived 
from  a stress  function  i|>  by  means  of 


x(l  - x) 

f(3x  - 2) (2  - x) 


x=a 


x— 3 


Fig.  14-.2 

Inequalities  1^.10. 
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X 


Fig.  14.3 

Lower  bounds  for  square  slab  with  square  cutout. 

uniaxial  tension 

biaxial  tension 
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dx=  9y2 


=4 

dx 


Q 

T £ w % 

xy  dxdy 


(15.1) 


then  the  equilibrium  equations  0.2  are  satisfied  identically. 
For  the  elastic  problem,  t|»  must  satisfy  the  biharmonic  equation 


= L^. 

3x 


+ 2^ 


. qS  - o 

p — p + — r u> 

ax  ay  ay 


(15.2) 


and  the  boundary  conditions.  However,  Eq.  15.2  is  derived  from 
Hooke's  law  and  does  not  apply  when  we  seek  only  a statically 
admissible  stress  field.  For  such  a purpose  we  may  choose  any 
function  r|>  which  satisfies  the  boundary  conditions. 

Once  is  chosen,  the  stress  components  are  determined 
from  15.1.  The  principal  stresses  are  then  given  by  [19,  Sec.  9 3 


(15.3) 


By  definitely  identifying  <5^  and  as  indicated  by  the  signs  in 
15.3,  the  yield  condition  1.1  then  states  that 


maxCd-p  - — 2k*  (15A) 

If  15. 4 is  everywhere  satisfied,  then  the  constructed  stress 
field  is  statically  admissible  and  will  furnish  a lower  bound  on 
the  load.  More  generally,  if  the  load  is  expressed  in  terms  of 
the  parameter  X,  any  value  of  \ which  satisfies  15* ^ is  a lower 
bound. 

The  stress  function  may  be  used  advantageously  with 
other  coordinate  systems.  Thus,  in  polar  coordinates  r,  0 the 
stress  components  defined  by 
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= 1 


d\j> 


6 = 


ar 

a2^ 

A 


+ -L  .0  , 
,2  80  » 


r© 


(15.5) 


will  satisfy  the  equilibrium  equations  identically.  The  princi- 
pal stress  are  then  given  by  equations  analgous  to  15. 4: 


cj 

1 


n 


+ (15.6) 


As  an  example,  let  us  consider  the  annular  slab  (Fig. 
12.1)  for  which  an  upper  bound  was  found  in  Sec.  12.  For  uniaxial 
tension  the  boundary  conditions  on  the  stresses  are 


at  r = a;  cJp  = 0,  = 

at  r = 1:  cJp  = k\(l  - cos 

In  view  of  the  dependence  of  the 
are  led  to  try  a stress  function 


(15.7) 

2©),  TrQ  = k\  sin  2©. 

boundary  conditions  upon  ©,  we 
of  the  type 


= f(r)  + g(r)cos  2© 


(15.8) 


where  f and  g are  to  be  determined.  It  then  follows  from  15*5 
that  the  stress  components  are 


_ f_*  + rg*  - 4g 
r ” r 


d = a—  + ~ — cos  2©, 
r2 


6 = f"  + g"  cos  2©, 

© 


(15.9) 


= 2 sin  2©, 


r© 


where  primes  indicate  differentiation  with  respect  to  r. 
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Substituting  15*9  into  15*7  we  obtain  the  boundary  conditions  for 
f and  g: 

f »(a)  = g(a)  = g*(a)  =0, 

(15.10a) 

f '(1)  = kX,  g(l)  = | kX,  g'(l)  = kX. 

Next,  let  us  consider  possible  discontinuous  stress 

fields  derived  from  stress  functions.  As  previously  discussed, 

the  normal  and  tangential  stress  must  be  continuous  across  any 

curve,  but  the  other  normal  stress  need  not  be.  For  the  annular 

slab,  15.8  implies  that  the  stress  function  is  analytic  as  a 

function  of  ©,  hence  the  only  possible  curve  of  discontinuity  is 

a circle  r = r,  . Across  such  a circle,  d and  t . must  be  con- 
k * r r© 

tinuous,  but  dQ  may  exhibit  a jump.  It  follows  from  15.9  that 
f and  g must  be  continuously  differentiable  but  may  have  discon- 
tinuous second  derivatives.  Thus, 

f'(r£)  = f'(r~),  g(r*)  = g(r“),  g'(r*)  = g'(r“).  (15.11) 

Subject  to  Eqs.  15.10  and  15.11,  we  can  choose  any 
functions  f and  g.  Once  these  are  chosen,  the  principal  stresses 
are  determined  by  15*9  and  15.6  and  X is  then  taken  to  be  the 
largest  number  such  that  the  yield  inequality  15. is  valid  at 
all  points  of  the  slab. 

The  difficulty  with  using  the  method  described  above  is 
that  even  when  the  functions  f and  g are  chosen  as  simply  as 
possible  so  as  to  satisfy  the  boundary  conditions  15.10,  the  re- 
sulting values  of  and  o2  become  quite  complicated,  so  that  an 
analysis  in  closed  form  becomes  impossible.  In  the  remaining  sec- 
tions of  this  chapter  we  shall  Indicate  two  possible  methods  of 
alleviating  this  situation  which  are  currently  being  investigated* 
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, 16.1 

16.  Numerical  techniques.  As  indicated  in  the  previous  sec- 

tion, even  a simple  choice  of  the  functions  f and  g in  Eq.  15.8 
will  lead  to  a complicated  function  of  r and  9 when  substituted 
into  the  yield  condition  15*4.  To  obtain  a lower  bound  for  the 
load,  it  is  necessary  to  find  that  point  in  the  slab  where  the 
left  hand  side  of  15*4  achieves  its  maximum.  However,  for  complex 
functions  this  cannot  be  done  in  closed  form  and  it  becomes 
necessary  to  use  numerical  techniques. 

A primary  disadvantage  of  numerical  methods  is  that  a 
separate  problem  must  be  solved  for  each  value  of  the  cutout 
radius  a.  However,  if  an  efficient  means  of  utilizing  modern  high- 
speed computing  machines  csn  be  devised,  this  drawback  becomes  less 
serious. 

In  the  present  section  we  shall  discuss  some  preliminary 
steps  which  have  been  made  in  this  connection,  and  indicate 
possible  extensions.  The  problem  has  only  recently  been  looked 
into,  and  the  results  and  suggestions  are  purely  tentative.  In- 
deed, it  is  impossible  to  say  whether  such  techniques  will  ever 
prove  useful. 

As  a first  step,  let  us  replace  f and  g in  Eq.  15.8  by 
two  new  functions  each,  defined  by 


f*(r)  = kXA(r)[l  + (1  - r)B(r)  ], 

2 

g (r)  = kXC(r)[l  + (1  -r)'D(r)]. 


(16.1) 


Substituting  16. 1 into  the  boundary  conditions  15*10,  we  see  that 
* The  results  of  this  section  were  obtained  by  E.  Levin. 
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B and  D are  entirely  arbitrary  (provided,  of  course,  that  they 
remain  finite  throughout  the  slab),  and  that  A and  C must  be  such 


that 

A(a)  = C(a)  * C«(a)  = 0, 

A(l)  = 1,  C(l)  = C*(l)  = 1. 


(16.2) 


For  any  A and  C satisfying  16.2,  we  express  successively 
the  stress  function  tjj(Bq.  15.8),  the  stress  components  dr,  oQt 
xrQ  ^s#  15.5),  the  principal  stresses  <s2  (Eq.  15*6),  and 
the  yield  condition  (l?*^)  in  terms  of  the  four  functions  A,  B, 

C,  D and  the  angle  0.  The  result  will  be  three  inequalities  which 
may  be  written 

XF^A,  B,  C,  D}  ©]  < 2, 


XF2[A,  B,  C,  Dj  ©]  < 2, 


(16.3) 


AF3[A,  B,  C,  D;  ©]  < 2. 

The  constant  X must  be  such  that  all  three  Inequalities  16.3  are 
satisfied  at  each  point  of  the  slab.  Any  choice  of  the  functions 
A,  B,  G,  D,  together  with  sufficient  computations,  will  furnish 
such  a X.  The  problem  then  is  to  choose  A,  B,  C,  D so  as  to  ob- 
tain the  largest  possible  value  of  X, 

Without  going  into  the  details  of  computation,  it  is 
evident  that  since  Cq  depends  upon  the  second  derivative  of  f and 
g,  the  values  of  F^,  Fg,  and  F^  will  depend  upon  B and  D,  even  on 
the  boundaries  r = a,  r = 1.  However,  a minor  change  in  the  defi- 
nitions 16. 1 will  eliminate  this.  Let  us,  then,  replace  16. 1 by 
the  definitions 
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f*(r)  = kXA(r)[l  + (1  - r)2(r  - a)B(p)], 

(16.4) 

g (r)  = kXC(r)[l  + (1  - r)3(r  - a)D(r)]. 

The  advantage  of  this  is  if  A and  C can  once  be  chosen  to  that  the 
maximum  value  of  2kP^  on  the  boundaries  appears  to  give  a reason- 
able value  for  X,  an  attempt  can  be  made  to  adjust  low  values  of 
in  the  interior  by  manipulating  B and  D.  Therefore,  in  decid- 
ing upon  the  '’best"  functions  A and  C,  it  is  necessary  to  consider 
only  values  at  the  boundaries  r = a,  r = 1. 

At  the  inner  boundary  r = a,  it  follows  from  the  bound- 
ary conditions  that  or  and  t are  both  zero,  so  that  the  yield 
inequality  becomes 

|oJ  <;  2k. 

Replacing  at  r = a by  its  value  in  terms  of  A and  C,  we  obtain 

- 2 < X(A*  + C‘»  cos  2©)  < 2.  (16.5) 


Inequality  16. 5 represents  a necessary  restriction  on  A and  C. 
Since  the  middle  member  is  linear  in  cos  20,  the  extreme  values 
are  at  © = 0 and  © = x/2.  If  the  assumed  stress  field  is  to  be 


at  all  close  to  the  actual  one,  dg  will  be  algebraically  greater 
at  © = 0 than  at  © = x/2,  so  that  C should  be  positive  at  r = a. 


As  an  example,  let  us  consider  the  particular  case 

a « -i,  and  take  simple  polynomials  satisfying  the  boundary  con- 

ditions:  , 

A = \ (3r  - 1),  B = 0,  D = 0, 


C = (3r  - ir(2  - r)/8. 


Since  only  a preliminary  analysis  is  Intended,  the  resulting 


Bll-22 


8>+ 


functions  Fk  are  evaluated  numerically  for  a coarse  mesh  in  r and 
©.  The  results  on  the  boundaries  are  given  in  Table  16.1 


© 

0 

30° 

450 

60° 

90° 

r = a 

Max  F^ 

5.25 

3.38 

1.50 

O.38 

2.25 

r = 1 

pi 

0.75 

1. 81+ 

2.06 

1. 77 

0.25 

r2 

0. 75 

1.73 

2.28 

2.57 

2.25 

*1 

0 

0. 11 

-0.22 

-0. 80 

-2.00 

Table  16.1 

Values  of  max  F^  on  the  boundary  for  Eq.  16.6. 

Since  the  maximum  value  of  F occurs  at  the  inner  boundary  at 
© = 0,  it  follows  from  16.3  that  a necessary  restriction  on  \ is 

^ - 0.  38.  ( 16#  7 ) 

That  tli is  stress  field  cannot  furnish  a "good”  lower 
bound  can  be  seen  as  follows.  Actual  plastic  yielding  cannot 
take  place  until  a plastic  region  extends  continuously  from  at 
least  one  point  on  the  inner  boundary  to  at  least  one  point  on 
the  outer  boundary,  since  otherwise  there  will  be  a determinate 
elastic  structure  remaining  which  can  sustain  a greater  load.  A 
necessary  condition  that  this  be  the  case  is  that  max  F^  has  the 
same  value  on  both  the  inner  and  outer  boundaries,  which  is  far 
from  the  case  in  Table  16.1. 

Therefore,  as  an  Improvement,  let  us  keep  the  same 
value  of  C,  but  take  A Identically  aero  on  the  inner  boundary. 
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Thus,  let  P 

A = (3r  - 1)  /4,  B = 0,  D = 0, 

C = (3r  - l)2(2r  - l)/8. 

The  values  on  the  boundaries  are  then  given  in  Table 


(16.8) 


16.2. 


e 

0° 

30° 

45° 

60° 

90° 

max  Fi 

3.75 

1.  88 

0 

1. 88 

3.75 

pi 

2.25 

2.74 

2,83 

2. 55 

1.75 

f2 

2.25 

2.93 

3.41 

3. 71 

3.75 

0 

-0. 19 

-0.59 

-1.16 

-2,00 

Table  16.2 

Values  of  max  F.^  on  the  boundary  for  Eq.  16, 8. 


In  this  case,  the  maximum  boundary  value  of  3.75  occurs  on  the 
Inner  boundary  at  0 = 0°  and  0 = 90°,  and  occurs  on  the  outer 
boundary  at  0 = 90°.  The  corresponding  necessary  restriction  on 
X is 

X £ 

In  Sec.  12  an  upper  bound  for  this  problem  was  obtained.  Setting 
a « Jj  in  Eq.  12,17  (or  reading  from  Fig.  12.2)  we  find  that 
X = 0.59  Is  an  upper  bound.  Since  16.9  is  reasonably  close  to 
that  upper  bound  it  would  appear  that  we  have  obtained  a reason- 
ably close  solution  on  the  boundary. 

However,  16.9  is  only  an  upper  bound  on  the  lower  bound. 
To  obtain  a true  lower  bound  for  the  assumed  function,  we  must 
compute  max  F*  throughout  the  slab.  The  results  are  shown  in 
Table  16.3. 


3^5 


= 0. 53. 


(16.9) 
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Nr\ 

9 

0° 

30° 

45° 

60° 

90° 

1/3 

“3.7? 

05 

d 

05 

OT" 

1/2 

3.38 

3.52 

3.56 

3.08 

1.88 

2/3 

3*30 

3.80 

3. 92 

3.26 

0,56 

5/6 

2.63 

4. 70 

5.11 

4.34 

1.88 

1 

2.25 

2.93 

3.41 

3.71 

3 .75 

Table  16.3 

Values  of  max  in  slab  for 

Eqs.  16.8 

• 

The  maximum  value 

for  the  mesh  considered  is  5.11,  occurring  at 

0 

= 45°.  Therefore, 

it  follows 

from  16.3 

that 

\ 

1 0.39. 

(16.10) 

The  next  logical  step  is  to 

keep  the  values  for  A and 

C given  by  Eqs.  16.8,  and  to  try  non-zero  values  for  B and  D in 
an  attempt  to  reduce  some  of  the  high  values  in  the  interior. 

Due  to  limitations  of  time,  we  are  unable  to  include  any  further 
results  of  this  technique  in  the  present  report. 

It  should  be  pointed  out  that  16.10  represents  only  a 
crude  approximation,  due  to  the  large  mesh  size  used.  However, 
it  appears  reasonable  to  use  a large  mesh  to  obtain  information 
as  to  '•good"  functions  A,  B,  C,  and  D.  Once  this  is  done,  the 
results  should  be  computed  for  a much  smaller  mesh  to  obtain  a 
more  precise  location  and  valuation  of  the  maximum  value  of  F^. 
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17.  Approximate  yield  conditions.  An  alternative  to  a completely 
numerical  approach  to  a statically  admissible  stress  field  is  to 
consider  an  approximation  to  the  yield  condition  sufficiently 
simple  so  that  analytic  techniques  may  be  employed.  It  follows 
from  Eqs.  1^.6  and  1|?#4  that  the  actual  yield  condition  consists 
of  the  three  inequalities 


u + 
- u + 


< 2k, 

< 2k, 

< k* 


(17. la) 
(17. lb) 
(17.  lc) 


where  we  have  introduced  the  notations 


u = (cp  + dg)/2,  v s (dp  - og)/2 . (17.2) 

If  we  sketch  the  surfaces  represented  by  the  equality  signs  in 
17.1  in  a Cartesian  space  with  coordinates  u,  v,  and  t,  we  obtain 
a surface  which  may  be  described  as  follows.  The  u axis  is  an 
axis  of  rotation;  for  |u|  < k the  surface  is  a cylinder  of  radius 
k;  for  k ^ |u|  £ 2k  the  surface  is  a pair  of  right  circular  cones 
capping  either  end  of  the  cylinder  with  vertices  at  u = + 2k. 

The  first  octant  of  such  a surface  is  sketched  in  Pig.  17.1. 

Since  we  are  concerned  with  finding  a lower  bound  on 
the  yield  load,  we  shall  consider  only  approximations  to  Fig.  17.1 
which  lie  wholly  within  the  actual  yield  surface.  If  this  is  done, 
then  any  lower  bound  obtained  using  an  approximate  yield  condition 
will  be  a true  lower  bound  for  the  problem# 
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T 


Yield  surface. 

The  difficulty  with  using  the  actual  yield  condition  is 
primarily  due  to  the  fact  that  the  stresses  contain  both  sin  2© 
and  cos  2©  (see  Eqs.  15*9)  and  hence  by  the  time  © is  eliminated, 
the  resulting  functions  of  r have  become  quite  complex.  There- 
fore, we  seek  an  approximation  which  will  simplify  the  dependence 
on  Q.  Since  is  the  only  stress  component  which  contains  sin 
2©,  while  u and  v contain  only  cos  20,  we  would  like  to  separate 
17*1  into  inequalities  containing  only  t and  those  not  containing 
*c.  Geometrically,  this  means  that  we  wish  to  Inscribe  a right 
cylinder  or  prism  with  generators  parallel  to  the  t axis’ in 
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Fig.  17.1. 

As  the  simplest  such  inscribed  prism,  consider  a hexag- 
onal prism  oriented  the  same  as  the  trace  of  the  true  yield  sur- 
face in  the  u,v  plane  and  bounded  by  the  planes  t = + k/\/2. 
Analytically  this  is  described  by  the  inequalities 

M < k/\/2  , jvj  < k/\/2, 

(17.3) 

The  choice  of  k/\/2  as  the  height  of  the  above  prism 
was  entirely  arbitrary  and  there  is  no  assurance  that  it  will 
produce  a reasonable  approximation.  To  obtain  a better  approxi- 
mation, let  us  consider  separately  the  cylindrical  and  conical 
portions  of  the  yield  surface.  In  the  cylindrical  portion  (Eq. 
17.1c),  it  is  sufficient  to  let  t and  v satisfy 

I v|  = k cos  p,  M < k sin  p.  ( 17. *+) 

In  the  most  general  case,  p can  be  a function  of  r and  Q,  but 
such  generality  would  destroy  the  simplicity  desired.  A reason- 
able basis  for  choice  would  be  to  divide  the  slab  into  a finite 
number  of  regions.  In  each  such  region  let  the  numerically 
largest  values  of  v and  v be  denoted  by  k\vn  and  k\xn  respective- 
ly. In  each  region,  P is  then  chosen  so  that  the  two  Inequalities 
17*4  are  equally  restrictive.  Thus,  for  the  correct  choice  of  P 

kXvn  = k cos  P,  k\Tn  = k sin  P, 
hence  p is  given  by 

P = tan-1  i , (17.5a) 

Tn 


and  the  cut-out  factor  must  satisfy 


X < sin  PAn.  (17.5b) 

In  the  conical  section,  we  construct  an  approximating 
prism  by  first  inscribing  a pyramid  in  the  cone,  and  then  a prism 
in  the  pyramid.  Since,  in  the  only  examples  considered  thus  far 
the  most  restrictive  part  of  the  yield  condition  corresponds  to 
the  cylinder,  the  weaker  approximation  under  the  cone  appears 
reasonable.  Analytically,  then,  we  note  that  17.1a,  b will  both 
be  satisfied  if 

|u|  + jvj  + H < 2k. 

Therefore,  we  demand  that 

ju|  + | v J < 2ka,  M £ 2k(l  - a).  (17.6) 


The  parameter  a is  determined  in  the  same  way  as  the  parameter  p 
in  Eq,  17.^.  In  particular,  if  a is  constant  throughout  some 


region,  then 
and 


a = 


% 


+ v 


n 


"n  + vn  + *n 


(17.7a) 


X <_22 . (17.7b) 

“ u_  + v_ 
n n 

If  Eqs.  15.9  are  substituted  into  17. 3 , each  of  the 
resulting  inequalities  are  linear  in  either  sin  29  or  cos  29, 
Since  all  inequalities  must  be  satisfied  for  all  values  of  0,  it 
is  thus  necessary  and  sufficient  to  satisfy  them  at  the  extreme 
values.  Therefore,  the  quantities  un,  vn,  xn  appearing  in  Eqs. 
17. 5 and  17. 7 are  independent  of  9, 

Let  us  assume  for  simplicity  that  each  of  the  regions 
is  a complete  annulus.  In  terms  of  the  functions  f and  g 
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The  possibility  of  investigating  the  annular  slab  by 
these  approximate  yield  conditions  is  only  now  being  investigated. 
Unfortunately  time  has  not  permitted  the  inclusion  of  any  examples 
of  the  technique  illustrated  by  Eqs.  17,  through  17,8,  In  the 
only  example  worked  out  so  far  by  means  of  Eqs.  17*3  (see  Sec.  21), 
the  resulting  lower  bound  is  disappointingly  low.  However,  at 
least  part  of  the  reason  appears  to  be  due  to  the  choice  of  stress 
function,  rather  than  to  the  use  of  an  approximate  yield  condi- 
tion. Further  investigations  of  this  approach  are  being  under- 
taken and  will  be  reported  on  elsewhere. 
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V.  EXAMPLES 

18,  Square  slab  with  square  cutout.  We  have  already  considered 
this  problem  in  uniaxial  tension  in  Secs,  13  and  14,  Therefore, 
it  follows  from  Sec,  3 that  we  need  only  find  a lower  bound  in 
equal  biaxial  tension  in  order  to  find  a bound  under  any  loading. 
As  in  the  case  of  uniaxial  tension,  it  is  useful  to  consider  a 
piecewise  constant  stress  field,  as  shown  in  Pig,  18. 1,  Due  to 
symmetry,  it  is  here  necessary  only  to  consider  one-eighth  of  the 
slab.  Thus  we  seek  to  determine  the  stresses  in  regions  1,  2, 
and  3 of  Fig,  18,1,  Symmetry  at  once  determines  the  value  of  0 
in  each  region: 


I 


ft 


Fig.  18.1 


Square  slab  in  biaxial  tension. 
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ex  = ©3  = o,  ©2  = £ , (18.1) 

while  the  boundary  conditions  show  immediately  that 


= 2kX,  s^  = 0.  (18.2) 

Applying  Eqs.  14.3  to  each  of  the  two  lines  of  discon- 
tinuity and  using  18. 1 and  18. 2 we  obtain 

(2k)(1  - cos  2ct^2)  + X(1  + cos  a12)  = - Xg  sln  2a  , 

- COS  2a  ) = u - V sin  2a  , 

23  2 ^ 23  (18.3) 

p 

(x  - 7jl)sin  2(^2  = Xg  cos  2a12, 

- (I3) sin  2a23  = x2  cos  2«23. 


These  equations  are  easily  solved  for  the  four  unknown  quantities, 
the  results  being 

2\  2X  tan  a, 


*2  ” tan 


13 


23 

- tan  a^,-,  ’ 2 tan  a,0  - tan  a,  7 


‘23 


‘12 


23 


(18.4) 


h , x 1 - *2 tan  g2.3. . la  = x 1 ; tan  X 

2k  tana^Ctan  al2-tan  a2y  2k  tan  a2^(tan  a12-tan  a2y 


The  two  angles  appearing  in  18.4  are  not  independent,  but  may 
both  be  expressed  in  terms  of  the  cutout  half-side  and  the  para- 
meter x.  Thus,  from  Fig.  18. 1, 


tan  a12  = 1 - x, 


tan  a23  = - (x  - a)/a* 


(18. 5) 


Substituting  18.5  into  18.4  we  obtain 
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■ 2Xa  /s\  x • & 

*2  = TTT-gS’  2 " 2 (1  - Hi  ’ 


(18.6) 


r*,  1 - x + a 

2k = x 


r3  . 2a  - x 

2k = " x 7x  "TTjrmn  • 


The  value  of  X must  now  be  determined  so  that  all  of 
14.4  are  valid.  The  governing  inequality  in  regions  2 and  3 may 
depend  upon  the  values  of  x and  a,  so  that  on  first  analysis  we 
must  solve  the  following  five  inequalities: 


II  - v 1 - * + a / -l 
2k  “ X (1  - x)(l  - a)  ~ lf 

r?  x • 2a 

2k  X (x  - aXl  - "a5 

r-3  2a  - x 

2k = x - 1’ 


(18. 7a) 
(18.7b) 

(18.7c) 


x2  _ 2\a 

2k  = x(l-  a) 


< 1, 

< 1. 


(I8.7d) 

(I8.7e) 


It  is  convenient  to  distinguish  two  cases,  depending  upon  the 
sign  of  x - 2a.  If  x - 2a  is  positive,  it  is  evident  that  all 
of  18.7  will  be  valid  if  18.7a  is.  However,  the  largest  X allow- 
able by  18. 7a  occurs  for  the  smallest  allowable  x,  which  in  this 
case  is  x = 2a.  Since  this  case  can  also  be  considered  under  the 
hypothesis  that  x - 2a  is  non-positive,  it  suffices  to  examine 
18.7  for 

x - 2a  £ 0. 

Inequalities  18. 7b  and  18. ?d  are  then  obviously  valid 
if  the  others  are,  so  that  we  seek  the  largest  X such  that 
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(18. 8a) 

(18. 8c) 

Considering  the  right  hand  sides  of  18. 8 as  functions  of  x,  it 
is  evident  that  all  three  are  monotone  and  positive  in  the  inter- 
val a <,  x £ min(2a,  1)  under  consideration.  Since  the  right 
hand  sides  of  18. 8a  and  18. 8c  are  decreasing  functions  of  x, 
while  18.8b  is  increasing,  the  situation  may  be  discussed  quali- 
tatively as  in  the  case  of  uniaxial  tension,  Fig.  14.2.  The 
optimum  choice  for  x will  be  the  intersection  of  18. 8b  with 
either  18. 8a  or  18. 8c,  whichever  yields  the  lower  value  of 
After  some  computation  it  can  be  verified  that 

x = 1 [1  + 2a  - \/(l  - a)2  + a2  ] , 

X = ■■  1* --  cV(l  - a)2  + a2  - a], 

1 - a 7 

for  ^ ( y'S  - 1)  < a < lj  (18.9) 

x = ^ [l  + 3a  - \J (1  - a)2  + 8a2]  , 

X 8 ik  Cl  + 3a  - V(1  - a)2  + 8a2], 

for  0 < a < ijf  ( \/5  - 1). 

The  resulting  function  of  \ is  shown  as  the  dashed 
curve  in  Fig.  14.3*  Since  the  solid  curve,  representing  the 
lower  bound  in  uniaxial  tension  is  everywhere  below  this  curve, 
it  follows  that  the  lower  bound  obtained  in  Sec.  14  (Eqs.  14.11) 
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Is  safe  for  any  loading. 

In  Fig.  18.2  we  have  redrawn  this  lower  bound,  together 
with  the  upper  bound  determined  in  Sec.  13  (Eqs*  13.13).  For 
a < A,  a better  upper  bound  is  given  by  the  method  of  Sec.  11#. 
Setting  Ty  = 1,  a0  = a in  Eq.  11. 3 the  latter  becomes 


X < 1 • cl. 


(18.10) 


Therefore,  for  a given  value  of  a,  the  cutout  factor X must  lie 
in  the  narrow  region  between  the  two  groups  of  curves  in  Fig. 
18.2.  For  practical  purposes,  we  would  appear  to  have  solved 
this  example  to  a reasonable  degree  of  accuracy. 

19.  Square  slab  with  a slit  [53.  In  the  previous  example  we 
were  able  to  obtain  a lower  bound  very  close  to  the  upper  bound 
by  using  a piecewise  constant  stress  field.  In  the  case  of  a 
square  slab  with  a slit  we  can  do  even  better  than  that  and 
obtain  the  actual  yield  load. 

The  discontinuous  stress  field  for  uniaxial  tension  is 
shown  in  Fig.  19.1.  The  method  of  solution  is  similar  to  that 
used  in  Sec.  l*f  and  the  results  are  as  follows  (for  details  of 
the  computations  see  [?]): 


©1  =0,  s-j_  = 2kX,  r^  = 2kXa/(l  - x), 

tan  2©2  = 2a/[x  + a(l  - a)], 

«2  = “ a(l  “ a)]/[a  + x(l  - a)],  (19.1) 

x2  = " X {,  + a(l  - a)]2  + *ta2^1</2/[a+x(l-a)], 

~ 0 j r ^ ~ 0 , s^  = 2kX/(l  “ a) , 

©1^  = 0,  3^  = 0,  r^  = - 2kXa/x. 


B 11-2 2 


Upper  and  lower  bounds 
on  yield  load  of  square  slab 
with  square  cutout. 
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2k  X. 


Fig.  19.1 

Square  slab  with  a slit-uniaxial  tension. 

Substituting  19.1  into  14.4  and  determining  the  governing  in- 
equality in  each  region  we  obtain  finally 

\a  <:  1 - x, 

\V  [x  + a(l  - a)]2  + 4a2  <:  a + x(l  - a), 

(19.2) 

\ < 1 - a, 

\a  < x. 

It  is  readily  verified  that  if  we  choose 

, 2 
x * 1 - a + a , 


(19.3) 
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then  X = 1 - a (19*^) 

will  satisfy  all  of  19«2  so  that  1 - a is  a lower  bound  on  the 
yield  load. 

That  1 - a is  also  an  upper  bound  for  uniaxial  tension 
follows  from  the  analysis  of  Sec.  11,  provided  we  assume  the 
idealized  uslit"  to  be  equal  in  width  to  the  thickness  of  the 
slab. 

Since  the  slab  is  not  symmetric  with  respect  to  its 
diagonals,  we  must  consider  both  biaxial  tension  and  uniaxial 
tension  in  the  x direction  in  order  to  show  that  1 - a is  the 
cutout  factor.  The  latter  is  trivially  true  since  in  this  case 
the  slit  has  no  effect  and  the  slab  would  support  a load  \ = 1, 

The  case  of  biaxial  tension  can  be  handled  with  a piecewise  con- 
stant stress  field  as  in  the  previous  example  (see  [51).  However, 
it  is  interesting  to  construct  an  entirely  different  statically 
admissible  stress  field  as  follows. 

We  divide  the  slab  into  three  regions  as  indicated  in 
Fig.  19,2.  The  region  interior  to  a circle  of  radius  a (region 
1)  we  take  to  be  stress  free,  while  the  regions  exterior  to  a 
unit  circle  (regions  3)  are  assigned  the  equal  biaxial  stress 
field. 

<3X  = <Jy  = dr  = <JQ  = 2k\,  t = 0 (19.5) 

Region  2 adapts  itself  naturally  to  polar  coordinates,  and  we 
consider  the  stress  field 

<jp  = 2k\(l  - a/r)/(l  - a), 


<jg  = 2kX/(l  - a),  t ^ = 0. 


(19.6) 
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Fig«  19.2 

Square  slab  with  a slit-biaxial  tension. 

Observe  that  19.6  furnishes  continuous  tractions  with  the  stress 
fields  in  regions  1 and  3,  and  satisfies  the  equilibrium  equation 

Jl.(rO  - d = 0 (19.7) 

dr  r © 

in  region  2.  It  remains  only  to  satisfy  the  yield  condition  1.1. 
Since  the  principle  stresses  are  in  the  coordinate  directions 
and  satisfy 

09.8) 

this  means  that  we  require 


d9  = 2kX/(l  - a)  < 2k, 


(19.9) 
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or 

X < 1 - a (19.10) 

In  agreement  with  19* 

The  preceding  types  of  stress  and  velocity  fields  can 
also  be  applied  to  slits  in  rectangular  slabs#  If  the  direction 
perpendicular  to  the  slit  is  the  longer  dimension,  we  again  ob- 
tain X = 1 - a as  the  cutout  factor.  Indeed,  the  same  upper 
bound  applies  as  before,  while  a statically  admissible  stress 
field  may  be  constructed  as  indicated  in  Pig.  19*3.  Region  1 
is  constant  stress 

°x  = 0,  «y  = 2k\,  = 0 


q.  Fig.  19.3 

Rectangular  slab  with  a slit. 
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while  region  2 is  the  same  piecewise  constant  stress  field  con- 
sidered for  the  square  slab# 

If  the  shorter  side  is  perpendicular  to  the  slit,  the 
upper  and  lower  bounds  obtained  no  longer  coincide#  A few  re- 
sults are  known  [20],  but  they  do  not  provide  good  approximations 
and  will  not  be  reviewed  here. 

20#  Square  slab  with  a circular  cutout  [3l«  An  important 
corollary  to  the  first  Prager-Greenberg-Drucker  theorem  is  that 
the  yield  load  of  an  elastic-perfectly  plastic  body  cannot  be 
lowered  by  the  addition  of  material  adjacent  to  a stress  free 
boundary.  Indeed,  it  is  merely  necessary  to  assign  a zero  state 
of  stress  to  the  new  material  and  complete  the  stress  field  as  if 
it  were  not  there#  It  was  this  idea  which  we  used  in  finding  a 
biaxial  lower  bound  in  the  previous  section. 

The  example  under  consideration  can  be  regarded  as  a 
square  cutout,  with  additional  material  added  to  change  it  to  a 
circle,  as  indicated  in  Fig.  20*1#  Therefore,  it  follows 


Fig#  20.1 


Square  slab  with  circular  cutout# 
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immediately  that  the  lower  bounds  obtained  for  the  square  cutout 
in  Secs.  14-  and  18  are  also  lower  bounds  in  this  case. 

The  solution  of  Sec.  13  also  provides  an  upper  bound 
for  the  present  example.  To  be  sure,  if  we  take  the  regions  of 
bending  to  have  any  finite  width,  then  Eq.  13,8  is  no  longer 
valid.  However,  in  Sec.  13  we  purposefully  chose  one  of  the 
angles  to  passing  to  its  limiting  value  of  zero,  and  left  the 
other  three  angles  ({3,  y»  anc*  & in  Fig.  13.2)  undetermined.  In 
the  present  application  let  all  four  angles  be  arbitrarily  close 
to  zero.  In  the  limit  we  again  obtain  Eq.  13. 11  for  the  inter- 
nal rate  of  dissipation  of  energy.  Since  the  external  energy 
dissipation  rate  is  unchanged  it  follows  that  Eqs.  13.13  are 
valid  in  this  case  also.  Therefore,  upper  and  lower  bounds  on 
the  yield  load  are  again  given  by  Fig.  18.2. 

21.  Annular  slab.  In  the  examples  previously  considered  we  have 
been  able  to  obtain  upper  and  lower  bounds  which  were  reasonably 
close  together.  However,  for  an  annular  slab  this  is  yet  to  be 
done,  and  the  best  upper  and  lower  bounds  obtained  so  far  differ 
considerably  from  each  other.  This  is  unfortunate  from  a practi- 
cal viewpoint,  because  probably  the  most  common  type  of  reinforced 
cutout  is  a circular  hole  with  an  annular  reinforcement.  There- 
fore, if  we  could  solve  the  annular  slab  problem,  we  could  use 
the  method  of  Sec.  4 to  solve  the  slab-independent  reinforcement 
problem.  Further,  as  we  shall  see  in  the  next  chapter,  the  only 
direct  methods  evolved  so  far  for  the  reinforced  slab  depend 
heavily  on  the  results  of  the  unreinforced  cutout  problem. 

In  view  of  the  importance  of  the  annular  slab,  it  seems 
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worthwhile  to  discuss  such  meager  results  as  have  been  obtained 
so  far,  inconclusive  as  these  may  be.  The  methods  used  have  been 
discussed  in  some  detail  in  previous  sections,  and  we  shall  here 
merely  collect  the  results* 

We  first  of  all  observe  that  the  case  of  equal  biaxial 
tension  can  be  easily  disposed  of,  by  considering  the  stress 
field  [1] 

<3  = 2kX(l  - 3)/(l  - a), 

r **  (21.1) 

<JQ  a 2kX/(l  - a),  TTrQ  a 0. 


Observe  that  this  is  the  same  stress  field  used  in  region  2 of 
Fig.  19.2  for  the  slab  with  a slit.  Therefore,  just  as  in  Sec.  19 
we  obtain  a lower  bound  of  X = 1 - a for  equal  biaxial  tension. 
Also,  using  sliding  out  of  the  plane  (Sec.  11)  we  obtain 
X = 1 - a as  an  upper  bound  so  that  the  cutout  factor  in  equal 
biaxial  tension  is 

X = 1 - a.  (21.2) 


It  remains,  therefore,  to  consider  uniaxial  tension. 
Upper  bounds  for  uniaxial  tension  can  be  found  by  the  methods  of 
Secs.  12  and  13.  Equations  12.15  and  12.17  provide  an  upper 
bound  which  is  less  than  1 - a for  all  values  of  a,  thus  showing 
that  the  limit  in  uniaxial  tension  is  definitely  less  than  in 
equal  biaxial  tension.  However,  for  large  values  of  a,  the  con- 
cept of  bending,  as  illustrated  in  Sec.  13  provides  a still  lower 
upper  bound.  Since  the  results  obtained  for  the  hollow  square 
slab  (Fig.  13.2)  were  independent  of  the  angles  a,  {5,  6 we 

may,  as  in  the  previous  section,  set  each  of  these  angles  equal 
to  zero  and  hence  Eqs.  13*13  are  valid  for  the  annular  slab  as 
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well.  As  a matter  of  fact,  since  for  the  annular  slab  an  in- 
crease of  a or  y from  zero  would  decrease  the  total  internal  rate 
of  dissipation  of  energy  in  regions  1 and  4,  (Fig.  13.2)  respec- 
tively, a still  lower  upper  bound  could  be  obtained  by  variation 
of  these  angles.  This  has  not  been  done,  however.  In  Fig,  21.1, 
the  upper  curve  shows  the  resultant  upper  bound  obtained  by  tak- 
ing the  lower  value  of  X as  given  by  Eqs.  12,17  or  13.13, 

For  any  stress  function,  a lower  bound  may  be  obtained 
using  Eqs.  17.3*  R»  K.  Froyd  has  used  these  equations  with 
stress  functions  defined  by 

f » = 0 g = Ci  + a)»  (r  - a)2,  a < r < 1-t 3 

2 (1  - ar  ~ “ 2 

f ' = [r  + "HP  1 (21.3) 

< r < 1 

g = M [l_n_a§  r2  + 2ar  - a ] 

2 (i  + ar  1 + a J 


Considering  only  the  case  a > ^ he  obtained  a lower  bound 


X < — , 

(1  + 3a) 


^ < a < \ , 


UlLtilH.:.42,  1 < a < jy-a-  --  S ± ^97-68/2  ^0>?8 
V2  a(3  - a)  2 bV2 

(21,  If) 

X.  < it, -=  1 6 \fl  - / 97-68  \/2  < a < 1. 

“2  l + a 4y^ 


The  resulting  lower  bound  is  sketched  in  Fig,  21.1, 

For  the  particular  case  a = i E,  Levin  has  considered 

3 

numerically  the  stress  function  defined  by 
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0 0.2  0.4  0.6  0.8  1.0 


Fig.  21.1 

Upper  and  lower  bounds  for  annular  slab. 
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f » 


g 


k\  (r  - -i) , 

3 

(k|)(3r  - l)2, 


f ' = kX  (r  - 1), 

3 

g = kX  (r  - i), 

f » = k\(3r  - 2), 

g = k X (r  - -i). 


J 


(21*5) 


The  worst  point  was  found  to  occur  at  0 = it  A,  r = 2/3*  The 
result  is  subject  to  some  question,  of  course,  as  a rather  course 
mesh  was  used.  However,  since  the  functions  involved  are  simple 
ones,  it  appears  reasonable  that  the  worst  point  would  be  at  the 
end  point  of  an  interval  in  r,  so  that  Levin's  result  may  be  ten- 
tatively accepted.  It  would,  however,  be  desirable  to  check  it 
analytically,  at  least  using  the  approximate  yield  condition. 

Levin's  result  is  considerably  higher  than  any  bound  so 
far  obtained  for  an  arbitrary  value  of  a.  Therefore,  for  a a 
little  less  than  1/3,  the  best  bound  available  at  present  is  to 
combine  21.5  with  a stress-free  region  for  r <1/3*  This  accounts 
for  the  horizontal  portion  of  the  lower  curve  in  Fig.  21.1* 

Finally,  for  a less  than  about  0.265,  the  best  lower 
bound  has  been  obtained  by  inscribing  a hollow  square  in  the 
annulus.  That  part  of  the  annulus  exterior  to  the  square  is  taken 
to  have  a constant  stress  state,  while  that  part  interior  to  the 
cutout  in  the  square  is  taken  to  be  stress-free.  Within  the 
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hollow  square,  the  discontinuous  stress  field  of  Sec.  14  is  used. 
Since  the  exterior  side  of  the  square  is  \/2  in  length,  while  a 
length  2 was  assumed  in  obtaining  Eqs.  14.11,  the  result  for  the 
present  case  is 

\A(a)  = X g ( a \/2 ) (21.6) 

where  X is  the  cutout  factor  for  the  square  as  given  by  Eqs. 
s 

14.11  (Eq.  14,  lib  is  the  applicable  one  in  the  range  of  interest 
for  a),  and  XA  is  the  desired  cutout  factor  for  the  annulus. 

Figure  21.1  shows  the  best  upper  and  lower  bound  cor- 
responding to  each  value  of  a.  As  previously  admitted,  the 
solution  is  by  no  means  satisfactory  and  the  investigation  is 
being  continued  in  the  hope  of  obtaining  better  bounds. 
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VI.  REINFORCED  CUTOUTS 

22.  Reinforcement  to  full  strength.  In  Sec.  4,  a method  was 
indicated  by  which  a reinforcement  for  full  strength  could  be 
designed  on  the  basis  of  the  analysis  of  the  hub  alone.  In- 
tuitively, one  might  expect  that  such  an  analysis  would  always 
be  conservative,  and  that  a consideration  of  the  stresses  in  the 
base  slab  would  lead  to  substantially  reduced  hub  dimensions. 
However,  as  we  shall  see,  this  is  not  always  the  case. 

As  an  example,  let  us  consider  a square  slab  of  half- 
side d with  a circular  cutout  of  radius  a,  reinforced  by  concen- 
tric circular  cylindrical  rings  of  outer  radius  b.  The  thickness 
of  the  hub  is  H,  and  the  slab  is  subjected  to  equal  biaxial  ten- 
sions 2k\(Fig.  22.1).  iy 
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Consider  first  the  radially  symmetric  stress  field  [1] 


d = 2k(l  - -S),  <5.  = 2k.  a < r < b 

r r ’ © ~ ~ 

«r  = 2k  |l  - |[1  -|(X  - |)]],  o#  = 2k,  b<r<l,  (22.1) 
d = 2k\,  d.  = 2k\,  r > b. 

Observe  that  this  stress  field  satisfies  the  equilibrium  equation 
(Bq.  19*7)  in  each  of  the  three  regions,  and  the  boundary  condi- 
tion = 0 at  r = 0.  Further,  at  the  hub  boundary  r = b the 
traction  per  unit  length  is 

2kH(l  - 3) 
b 

computed  from  either  side.  Finally,  the  tractions  applied  at 
r = 1 will  be  in  equilibrium  provided  that 


X = 1 - i + (|)(b  - a)/d.  (22.2) 

a h 

Therefore,  Eqs.  22.1  furnish  an  equilibrium  stress  field.  If  we 
impose  the  additional  condition  that 


(22.3) 


then  at  all  points  of  the  reinforced  slab 


°v  1°  (22.4) 

so  that  the  stress  field  22,1  is  statically  admissible.  It 
follows  from  Theorem  1 that 

X 2 1 - | + (|)(b  - a)/d,  (22.5) 

Next,  we  may  construct  a kinematically  admissible 


velocity  field  based  on  sliding  out  of  the  plane,  as  in  Sec.  11, 
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An  elementary  computation  then  shows  that  1 - b + (H/h)(b  - a) 
is  also  an  upper  bound  for  the  yield  load  X.  Therefore,  the  load 
and  hub  thickness  are  related  by  Eq,  22.2. 

If  the  reinforcement  is  to  be  designed  for  full  strength, 
then  we  may  set  X = 1 in  Eq.  22,2  and  obtain 


(22.6) 


This  result  is  independent  of  the  dimension  d which 
characterizes  the  slab.  Indeed,  the  same  result  could  have  been 
more  easily  obtained  by  the  method  of  Sec.  4.  For  a plane  annular 
slab  of  inner  and  outer  radii  a and  b respectively,  the  yield 
load  under  biaxial  tension  is  easily  shown  to  be  [1], 


X=  1-f  . (22.7) 

Substituting  this  into  4.1  yields  22.6  Immediately. 

Although  reinforcement  to  full  strength  is  independent 
of  d,  the  same  is  not  true  if  we  wish  to  reinforce  to  some  value 
less  than  full  strength.  Solving  22.  5 for  H/h,  we  obtain 


M - JLrJk 

h b - a b - a 


d. 


(22.8) 


In  Fig,  22.2  we  have  sketched  H/h  as  a function  of  d for  the 
particular  case  a = 1,  b = 2 for  different  values  of  X,  Observe 
that  for  any  value  of  X < 1 a sufficiently  large  value  of  d brings 
the  slab  to  the  required  strength  without  any  reinforcement 
(H/h  = 1). 
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H/h 


Fig.  22.2 

Reinforcement  to  less  than  full  strength. 

23 o Wide  square  reinforcement.  As  a second  example  of  direct 
computation  of  reinforced  cutouts,  consider  a square  reinforce- 
ment  of  a square  coutout  in  a square  slab  (Fig.  23.1).  The  state 
of  stress  in  the  hub,  as  well  as  in  the  base  slab,  Is  to  be 
approximated  by  plane  stress. 

Let  us  consider  a piecewise  constant  stress  field, 
dividing  the  base  slab  into  four  regions  as  in  the  previous  sec- 
tion, and  dividing  the  hub  into  four  similar  regions  (Fig.  23.1), 
In  each  of  regions  1>  3»  5»  7 and  8,  symmetry  demands  that  the 

principal  directions  be  the  coordinate  directions  so  that 
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i 

I 

Fig.  23.1 


Wide  square  reinforcement  of  square  slab 
with  square  cutout. 
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«1  = «3  = e4  = ©5  = ©7  = ©g  = o.  (23.1a) 

Equilibrium  with  the  applied  boundary  tractions  determines 


s 


1 


2k\, 


= 0 


> 


'8 


= 0. 


Next,  if  we  denote  d in  region  3 by  2ty.,  and  define 


n 


(23.1b) 


as  the  ratio  of  base  slab  height  to  hub  height,  equilibrium 
between  regions  3 and  5 furnishes 


s^  = 2ky,  = 2kp.q,  (23.1c) 

In  an  example  such  as  the  present  one  which  contains  several 
regions,  it  is  desirable  to  linearize  the  equilibrium  equations 
between  adjoining  regions  (Eqs.  14.3)  by  defining 


u = x cos  20,  v = x sin  20. 


(23.2a) 


Using  the  further  notation 

p = , o = -2- 

p 2k»  2k 

the  situation  in  each  region  is  given  in  Table  23,1. 

Since  we  have  already  satisfied  the  boundary  conditions 
at  the  inner  and  outer  boundaries  and  between  the  hub  and  base 
slab,  there  are  a total  of  twelve  Eqs.  14.3  as  follows; 
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Region 

© 

P 

U 

X 

u 

V 

1 

1 

0 

pi 

X 

p + X 
1 

pl- 

X 

pl-x 

0 

2 

©2 

p2 

d2 

W2 

X2 

u2 

v2 

3 

0 

p3 

P 

P3  + P 

p3  - 

P 

P3  - n 

0 

t 

4 

0 

0 

d^f 

d!+ 

"d4 

-°4 

0 

? 

0 

p5 

pn 

P5  + pn 

p5- 

pn 

p5  - pn 

0 

6 

P6 

d6 

w6 

*6 

u6 

v6 

7 

0 

p7 

0 

p7 

p7 

P7 

0 

8 

0 

0 

d8 

°8 

-a8 

-d8 

0 

Table  23. 1 

Stress  quantities  in  Pig.  23.1. 

regions 

1,2: 

(Pi  + 

X)  - 

(P]_  " ^cl  = 

w2  “ 

u2Cl 

- v2S-j_, 

(Pi  - 

Ms1 

r 

u2S1 

- v2Clf 

regions 

3,2: 

(p3  + 

fl ) - 

(p3  - M-  )C3  = 

W3  ' 

U2C2 

- V3> 

<p3  - 

p)s3 

- 

u.S 
3 3 

- v2C3, 

regions 

4-, 2: 

\ “ 

W2  " 

v2> 

-a4  = 

u2, 

regions 

5,6: 

Cp5  + 

pn)  - 

(p  5 - pnX5= 

w6  ' 

U6C5 

- V6S?» 

(p5- 

pn)s^ 

rr 

“6S5 

" V6C?> 

regions 

7,6: 

P7U 

- C7) 

= 

w6  " 

«6°7 

- V6S7> 

P7S7  = 

“6®7 

- v6C7, 
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regions  8,6:  °8  ~ ~ v6» 

**  d8  = u6» 

where  we  have  used  the  symbols 

= cos  2a12,  = sin  2a  , 

C-.  = cos  2a  , S~  = sin  2a  , 

3 32’  3 32J  (23.4) 

C5  = cos  2a?6,  = sin  2a?6, 

C?  = cos  2a76,  S?  = sin  2a^ 

Equations  23.3  separate  into  two  groups  of  six  equations 
each,  the  first  six  equations  containing  p-p  p^»  w2,  u2, 

and  V£  as  unknowns;  the  second  group  containing  p^,  p^,  dg, 

u^  and  v^.  The  equations  are  easily  solved  with  the  use  of 
the  trignome trie  identities 


cos  2a  = (1  - tan2a)/(l  + tan2a), 

2 

sin  2a  = 2 tana/(l  + tan  a). 
Further,  it  follows  from  Fig,  23,1  that 

y 

tan  a^2  = 1 - y,  tan  a^2  = 1 - 

tan  «?6  = 1 - |,  tan  a?6  = 1 - 

After  some  computation,  then,  we  finally  obtain 


(23.5) 


(23.6) 
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Ug  = 


X - |xb 
" 1 - b ’ 


v2  = 


P?  = 


P7  = 


-gX,-  u b 

y 1 - b 
ab  M-  n 

rr^  • tT^x 

^ 

b - a a - x 


°8  = 


b - a 


> 


u - 

6 


x - 2a  b 
x b - a 


u6  = 


b - a 


v*  = 


- la  b 
x b - a 


nn< 


(23.7) 


Next,  we  compute  the  governing  inequality  in  each 
region.  However,  we  first  record  the  elementary  inequalities 


0 < a<x<b<y<l, 

0 < rj  < 1,  (23.8a) 


0 < M-  < X < 1. 

Taking  23.8a  for  granted  wherever  they  are  needed,  we  see  that  in 
region  1,  p and  6 are  both  positive  and  0 < 1 is  included  in 
23.8a  so  that  the  governing  inequality  is 


= X » u 


< 1. 


'1  1 - y 1 - b — 

In  region  2,  | x I is  easily  computed  from  the  identity 


(23.8b) 


lx!  = U + v2, 

and  it  can  be  shown  that  for  any  values  of  the  parameters 
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I X I > |w|.  Therefore,  the  governing  inequality  is 

lx2l  \/(x  - v*?  + C^>Z(>-  - n)2  < 1.  (23.80) 

In  region  3,  on  the  other  hand,  we  have  0 > 0 > p for  any  values 
of  the  parameters,  so  that  the  governing  inequality  is 

!x3I  = a3  - P3  = H + y-Tb  T^b  - X‘  (23-8d) 

Since  p = 0 in  region  4,  the  governing  inequality  is 

d = iff  < 1.  (23,8e) 

4 1 - D — 


Finally,  similar  arguments  for  each  of  the  four  regions  of  the 
hub  lead  to  the  governing  inequalities 


* _ abug 

5 ~ (b  - a)(b 


(23. 8f) 

(23. 8g) 

(23. 8h) 

(23.81) 


Inequalities  23.8  may  be  considered  in  either  of  two 
ways.  On  the  one  hand,  for  given  values  of  a and  b we  may  wish 
to  determine  the  maximum  value  of  q(=  h/H)  which  will  restore 
the  slab  to  full  strength.  As  far  as  the  present  type  of  stress 
field  is  concerned,  the  situation  is  similar  to  that  encountered 
in  the  previous  section  in  that  we  obtain  no  Improvement  over  the 
result  obtained  by  applying  the  theory  of  Sec.  4 to  the  example 
considered  in  Sec,  18,  To  see  this,  we  need  only  consider  In- 
equality 23,8,  With  \ = 1,  this  becomes 
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1 - jib  < 1 - b. 

Since  b must  be  positive  and  p < 1 it  follows  that  p = 1.  There- 
fore, the  applied  boundary  tractions  are  transmitted  directly  to 
the  hub,  as  stated. 

A second  problem  is  to  find  the  yield  load  of  a given 
reinforced  slab.  That  is,  given  a,  b,  and  t),  find  the  maximum 
X satisfying  Inequalities  23.8.  Since  there  are  three  parameters 
in  the  problem,  it  is  impracticable  to  obtain  any  general  results. 
However,  a solution  is  easily  obtained  for  any  particular  value 
of  the  parameters. 

As  an  example,  consider  a reinforced  slab  with  a = ^ , 
b = q = -i  (observe  that  for  a circular  hole  this  would  provide 
full  strength).  Since  Inequalities  23. 8e  and  23. 8i  are  conse- 
quences of  23.8c  and  23. 8g,  respectively,  we  must  consider  the 
six  inequalities 


X - p < 1 - y, 

(23.9b) 

WX  - |)2  + (X  - p)2/y2]  < 1, 

(23. 9c) 

P + (X  - p)/(y  - -i)  <1, 

(23, 9d) 

P < 2 - 4x, 

(23. 9f) 

M-  \ A-  + -“^5  < 1» 

V 4x2 

(23. 9g) 

p <kx  - 1. 

(23. 9h) 

Since  p < X,  it  is  obvious  that  a larger  value  of  p 
will  always  allow  a larger  value  of  X in  Inequalities  23.9b,  c. 

An  elementary  calculation  shows  that  the  same  result  is  valid  for 
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23. 9d.  Since  these  are  the  only  three  inequalities  which  contain 
X and  y,  while  the  remaining  three  are  the  only  ones  which  con- 
tain x,  it  follows  that  we  wish  to  determine  x in  23. 9f,  g,  h so 
as  to  maximize  in  those  inequalities.  An  analysis  similar  to 
that  used  in  Sec#  lb-  shows  that  we  must  choose 


x 


„ 3 


(23.10) 


With  this  value  of  p.,  the  first  three  Inequalities 
23.9  become  ~ 

X<|-y,  (23.11b) 

2 \f\ T^)2  + (X  - l)2/y2  < 1,  (23.11c) 

X .1  y + • (23.  lid) 

These  inequalities  are  again  similar  to  those  considered  in  Sec. 

lb.  Let  us  try  the  value  of  y for  which  23.11b,  d are  both 

equalities.  Then  a 0 

y = X = (23.12) 

It  is  readily  verified  that  these  values  also  satisfy  23. lie  so 
that  we  have  obtained  the  best  value  of  X for  the  considered 
stress  field. 

By  way  of  comparison,  the  method  of  Sec.  lb-  together 
with  the  results  of  Sec.  lb-  would  furnish  a lower  bound  of  X = 

The  upper  bound  obtained  by  considering  sliding  out  of 
the  plane  is  easily  seen  to  be  X = 1,  corresponding  to  full 
strength.  However,  by  using  a bending  type  velocity  field  as  in 
Sec.  13,  we  can  show  that  full  strength  is  not  attained.  For 
simplicity,  we  shall  assign  convenient  values  to  the  hinge  centers- 


i 
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by  varying  these  values  we  could  presumably  reduce  the  upper 
bound  still  further.  Thus,  with  reference  to  Fig.  13.2,  take 
n o = o,  f Observe  that  this  locates  the  top 

hinge  at  the  junction  of  hub  and  base  slab,  and  the  lower  hinge 
at  the  outer  edge  of  the  slab.  The  internal  energy  dissipation 
rates  are  given  by  Eqs.  13.6,  where  h must  be  replaced  by  H = 2h 
in  areas  corresponding  to  the  hub.  The  resulting  total  is 

°®i  = [(tf}  + Cj)  + 0 + = 15  kh0/l6.  (23.13) 

Since  the  external  rate  of  dissipation  of  energy  is  still  given 
by  13«10»  we  have 

= 2k\h0(£  + ji  - £)  = khOX.  (23.14) 

Comparing  23. 13  and  23, 14  we  obtain  the  upper  bound 

(23.1?) 

24.  Wide  slabs  with  narrow  reinforcements.  If  the  hub  is  suffi- 
ciently narrow  to  be  approximated  by  a curved  beam,  as  in  Chap. 

II,  then  the  methods  of  Chaps.  II  and  IV  can  be  combined  to  fur- 
nish a lower  bound. 

As  an  example,  let  us  consider  a circular  cutout  with 
an  annular  cylindrical  reinforcement  subject  to  uniaxial  tension. 
We  divide  the  base  slab  into  several  regions  as  indicated  in 
Fig.  24.1.  Regions  1 through  5 constitute  the  base  slab  and 
region  6 is  the  hub.  The  stress  distribution  in  the  hub  is  ex- 
pressed in  terms  of  stress  resultants,  as  in  Chap.  II.  The  stress 
distribution  in  regions  1 through  4 is  similar  to  that  discussed 
in  the  base  slab  in  the  previous  section.  The  component  of  stress 
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d is  2\k  in  region  1 to  be  in  equilibrium  with  the  applied  load, 

y 

while  in  region  3 it  is  2k p where  (i  is  a constant  to  be  deter- 
mined. In  region  5,  the  stress  is  determined  by  regions  3 and  4 
to  be 

dx  = °>  °y  = 2kPf  Txy  = °* 

Therefore,  it  follows  that  the  hub  is  acted  upon  by  a uniformly 
distributed  vertical  force  of  magnitude  2kph  per  unit  length. 

The  forces  acting  upon  regions  1 through  4 are  exactly 
the  same  as  those  acting  on  the  corresponding  regions  in  Pig. 
23.1,  considered  in  the  previous  section.  Therefore,  the  same 
stress  distribution  will  subsist  and  hence  the  same  governing 
inequalities.  From  Eqs.  23*  8b  through  e we  have  then 


2LZJ*  b <± 

1 - y 1 - b - 9 

(24. la) 

(\  - pb)  + (^)  (X  - p)^  <;  (1  - b)^, 
¥ 

(24. lb) 

“ + rb*1’ 

(24.  lc) 

X - pb  < 1 - b. 

The  last  of  these  inequalities  is  a consequence  of  24,1b  and 
hence  will  not  be  considered  further.  Region  of  course,  does 
not  add  any  new  inequalities*  The  inequalities  in  region  6 are 
obtained  from  Eqs.  6.12,  replacing  \ by  p*  Thus 

rh*  <*-*> 

4<$  . (24.2b) 

Next,  we  observe  that  the  left  hand  sides  of  24.1a  and 
c are  decreasing  functions  of  p for  all  values  of  y and  b,  while 
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the  same  will  be  true  for  24* lb  if 

y < b(2  - b)/(l  - b).  (24.3) 


It  is  found  that  this  is  generally  the  case,  so  that  we  shall 
assume  that  the  left  hand  sides  of  all  three  Inequalities  24.1 
are  decreasing  functions  of  +l«  It  follows  that  since  p is  bound- 
ed only  from  above  by  inequalities  24.2,  we  should  choose  p to  be 
the  largest  value  permitted  by  these  equations.  Then  p becomes 
a known  quantity  and  inequalities  24.1  are  to  be  solved,  deter- 
mining y so  as  to  maximize  X,  It  should  be  pointed  out  that  if 
b and  y should  turn  out  to  be  such  that  24,3  is  reversed,  this 
will  not  lead  to  an  invalid  result, since  for  anv  choice  of  p and 
y a legitimate  lower  bound  on  X is  obtained.  It  merely  means 
that  the  best  possible  lower  bound  for  the  given  data  would  not 
be  obtained* 

With  p fixed,  then,  the  general  character  of  inequal- 
ities 24,1  is  the  same  as  that  illustrated  in  Fig.  14. 2 for  in- 
equalities 14,10.  Thus,  depending  upon  the  values  of  the  para- 
meters, the  best  value  for  y will  be  the  intersection  of  24.1c 
with  either  24,1a  or  24,1b.  After  some  computation  it  is  found 
that 

X > minCX^,  X2),  (24,4) 

where 

\ = V + , (24.5a) 


and  X2  is  defined  implicitly  by 

,2 


1 - 


(X2  - p)b' 
1 - b 


4b*(X2  - p)‘ 


(1  - b)2  - (X.  - bn)* 


(24.  ?b) 
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A convenient  upper  bound  for  the  above  example  is  ob- 
tained by  considering  sliding  out  of  the  plane: 


\<l-b  + (b-  a)  (|-). 


(24,6) 


Other  upper  bounds  can  also  be  obtained  by  considering  bending  as 
in  the  previous  section.  Some  specific  examples  will  be  given  in 
Sec.  26, 

The  same  approach  can  be  used  for  any  shaped  hub.  In 
each  case  p,  is  obtained  by  the  methods  of  Chap.  II,  and  X is 
then  given  by  24.4,  Thus,  for  a square  cutout  with  a square 
reinforcement,  it  follows  from  6. 18  that 


4 


min  -S 
h 


\/2  (a  + 6 ) 


H -a  + V a^  + 4&2 
h 2 (a  + bl 


(24.7) 


As  an  example,  let  us  consider  the  reinforced  slab  for  which  a 
plane  stress  bound  was  found  in  Sec.  23.  With  a = 6 = 1/4, 

H/h  = 2 Eq.  24 .7  furnished  p,  = v/j?  - 1)  - 0.62,  Substituting 
this  into  24.5a  we  obtain  X-^  = 0,76.  It  is  easily  shown  that  X2 
is  greater  than  X^,  in  this  case,  so  that 


X = 0.76 


(24.8) 


is  a lower  bound.  Comparing  this  with  Eqs.  23.12  and  23. 15  we 
see  that  this  type  of  approach  gives  a value  between  the  upper 
and  lower  bounds  furnished  by  the  plane  stress  approach. 


BL1-22 


12? 


25.  Other  loading  conditions.  Up  to  now  we  have  been  consider- 
ing only  the  situation  where  a uniform  load  has  been  applied  to 
the  boundary.  However,  in  applications,  it  may  be  that  only  the 
resultant  load  is  known,  and  not  the  load  distribution*  In  the 
present  section  we  shall  discuss  this  situation  more  generally. 

In  the  first  place,  we  observe  that  an  assumption  must 
be  made  that  the  load  intensity  per  unit  length  is  nowhere  greater 
than  2kh,  since  otherwise  the  material  of  the  slab  would  fail 
locally*  Secondly, it  is  evident  that  the  slab  is  generally  stronger 
if  the  load  is  applied  near  the  edge  of  the  slab,  away  from  the 
cutout.  Therefore,  we  may  consider  two  extreme  cases  as  indicated 
in  Fig.  25.  1.  In  each  case,  a vertical  load  of  intensity  2kh 
per  unit  length  is  applied  over  a length  c in  each  quadrant.  In 
Fig,  25.1a  the  load  is  at  the  outer  edge,  and  in  Fig,  25.1b  at 
the  center.  In  either  case  we  may  define  a design  problem:  to 
determine  a reinforcement  which  will  support  a given  load  span  c; 
or  an  analysis  problem:  to  determine  the  maximum  c for  a given 
slab. 

Let  us  consider  first  an  unreinforced  slab  with  an 
outer  load  applied.  For  the  case  c = 1 - aQ  a statically  admis- 
sible stress  field  may  be  trivially  constructed  as  indicated  in 
Fig.  25.2,  Region  1 is  stress  free,  while  in  region  2 


<3  = t = 0 

y xy 


, 3X  = 2k. 


Therefore,  it  follows  that  1 - aQ  is  a lower  bound  for  c.  On  the 
other  hand,  a kinematically  admissible  velocity  field  consisting 
of  sliding  out  of  the  plane  shows  that  1 - a0  is  also  an  upper 


bound,  so  that  in  this  case 


(b) 

Fig.  25.1 

Extreme  load  distributions. 
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Fig.  25.2 


Statically  admissible  stress  field  for  extreme  load. 

c = 1 - a0.  (25.1) 

As  an. application  of  this  type  of  analysis,  consider  a 
test  slab  which  is  clamped  perfectly  rigidly  at  either  end  (Fig. 
25.3).  If  the  clamp  is  perfectly  rigid  and  indefinitely  strong 
then  anv  equilibrium  stress  field  in  the  clamp  is  statically 
admissible.  Such  a stress  field  can  always  be  found  for  any  pre- 
scribed boundary  stresses,  for  example,  the  elastic  solution. 


| Fig.  25.3 

Statically  admissible  stress  field  for  clamped  slab. 
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Therefore,  in  Fig.  25.  3>  let  the  stress  field  in  the  clamp  be  such 
that  it  is  in  equilibrium  a total  load  of  2k\  on  AB  (in  whatever 
manner  it  is  actually  applied),  no  load  on  EC,  and  a uniform  load 
2k\/(l  - aQ)  on  CD.  The  other  sides  are  stress  free  and  there 
are  no  shearing  tractions.  Region  1 is  taken  as  stress  free, 
while  in  region  2 

dy  = ^xy  = °>  dx  = 2kK/(1  “ a0>« 

The  above  stress  field  will  be  statically  admissible  for  the 
entire  assembly  of  clamp  and  slab  provided  that  o in  region  2 is 
less  than  the  yeild  stress.  Thus  a lower  bound  for  \ is  1 - a0. 
Since  it  follows  from  Eq.  11. 3 that  this  is  also  an  upper  bound, 
we  have 

^ = 1 - a0.  (25.2) 

Therefore,  for  a perfectly  rigid  and  strong  clamp  the  collapse 
load  is  determined. 

In  actual  practice,  the  clamp  may  not  be  perfect,  so 
that  25.2  is  not  exact.  While  conceivably  the  actual  load  as  it 
reaches  the  slab  could  be  as  indicated  in  Fig,  25*  lb,  in  practice 
this  seems  most  unlikely,  so  that  the  true  collapse  load  is  prob- 
ably somewhere  in  between  25.2  and  the  uniform  load  analysis  con- 
sidered in  the  preceding  section. 

For  a reinforced  cutout  the  analysis  is  not  quite  so 
trivial.  As  a matter  of  fact,  for  reinforcement  to  full  strength 
the  present  viewpoint  offers  no  improvement.  However,  for  analy- 
sis under  a given  load  less  than  full  strength,  modifications  of 
the  above  procedure  can  certainly  be  used  to  advantage. 

As  an  example,  consider  the  reinforced  slab  discussed 
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in  the  previous  section,  pulled  in  uniaxial  tension  by  a perfect 
clamp  (Pig.  25.4).  The  stress  distribution  in  the  slab  is  the 
same  as  in  Fig.  25*3,  where  the  length  c is  chosen  so  that  the 
tractions  transmitted  across  AB  onto  the  slab  are  2kh  per  unit 
length.  The  base  slab  stresses  are  then  statically  admissible, 
and  it  remains  only  to  consider  the  hub. 

The  hub  is  subjected  to  a horizontal  force  of  magnitude 
2kh  per  unit  length  for  01  < © < %/2  and  to  no  force  for  0 < 

The  resultants  may  be  computed  by  equations  analogous  to  5*8, 
and  it  is  found  that 


N(Q)  = 0,  M(0)  = M(O),  0 <9  <QV 

H(9)  = 2fchb(  sin  9 - sin  9^  sin  9, 

M(©)  = M(0)  + khb(sin  © - sin  ©^)(a  sin  © 


(25. 3) 

< © < rc/2. 


Using  the  notations 


C = sin  ©,  £ = sin  ©^ 

r\  = h/H,  a = a/b,  (25.4) 

Y = 2HM(0)/kh2b2, 


and  substituting  25.3  into  6.3  we  obtain  the  yield  inequalities 

|I  + (Sx*;  - tjUaC  - Cx>l  + (C  - S (1  - o2)/n2»  (25.5) 

Inequality  25.5  must  be  satisfied  for  all  ^ rc/2,  and  Y is 

to  be  chosen  so  that  ^ is  a minimum.  It  follows  from  Pig.  25*4 
that  the  resultant  load  on  the  quarter  is  given  in  terms  of  by 
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2kXh  = 2kXc  = 2kh(  1 - b^), 


so  that  the  reinforced  cutout  factor  is 


X = 1 - bC 


1* 


(25.6) 


The  solution  of  inequalities  25.5  is  not  simple,  even 
numerically  for  specified  values  of  r\  and  a and  no  examples  have 
been  computed  as  yet.  Possibly  some  of  the  techniques  for  the 
limit  analysis  of  arches  [17]  can  be  used  to  advantage. 

A somewhat  simpler  (although  in  general  not  as  good) 
lower  bound  can  be  obtained  by  assuming  full  yield  stress  in 
region  2a  (Fig.  25. *f)  as  before,  but  a constant  traction  2khp.  in 
regions  1 and  2b.  The  hub  analysis  is  then  identical  with  that 
in  Sec.  6,  except  that  X is  to  be  replaced  by  p..  Thus,  from 


Eq  s . 6 * 12 , 


a = _H  6 

h(a  + 5) 


if  - > 2, 
a — 7 


p = H 


a..+  2t.  -.s  if  - < 2, 

h(a  + b)  a 


(25.7) 


Since  the  total  load  is 


2khX  = 2kh( 1 - a - 6)  + 2kh^(a  + 5), 


the  cutout  factor  is 


X = 1 - (a  + 6)  + — 

n 


if  - > 2, 

a — 


X = l - (a  + 6)  + [ V a2  + 2 b2  - a]/rj  if  - < 2. 

a ~ 


(25.8) 
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26.  Comparison  with  experiment.  The  only  known  experiments  on 
the  yield  loads  of  slabs  with  cutouts  were  carried  out  at  the 
University  of  Washington  by  Vasarhelyi  and  Hechtman.  These  ex- 
periments are  described  in  [21],  and  it  is  of  some  interest  to 
compare  the  experimentally  observed  results  with  the  theoretical 
results  here  obtained. 

A total  of  23  specimens  were  tested  involving  slabs 
without  cutouts,  with  cutouts,  and  with  reinforced  cutouts.  A 
variety  of  cutout  and  reinforcement  shapes  were  considered.  In 
each  case,  the  slab  was  36"  in  width,  and  had  a gauge  length  of 
36”  so  that  the  base  slab  was  essentially  square.  All  specimens 
were  lA"  thick,  and  were  made  of  a plain-carbon  semi-killed  grade 
steel  in  an  as  rolled  condition. 

We  shall  here  consider  only  those  slabs  with  circular 
cutouts  reinforced  with  concentric  annular  rings,  although  the 
methods  of  this  report  could  be  used  to  obtain  a theoretical 
result  for  all  the  specimens  tested.  In  order  to  avoid  possible 
experimental  difficulties  in  determining  the  yield  stress,  we 
shall  identify  the  cutout  factor  as  the  ratio  between  the  yield 
loads  of  the  slab  considered  to  the  same  base  slab  without  cutout 
or  reinforcement.  Further,  we  shall  define  the  yield  load  as  the 
load  for  which  a marked  increase  in  deformation  rate  first  occurred. 
In  [21]  this  is  called  "general  yielding".  In  each  case  the  slab 
continued  to  sustain  greater  loads,  indicating  the  existence  of 
strain-hardening.  Finally,  we  shall  use  the  hub-base  slab  analy- 
sis, regardless  of  the  actual  welding  process  by  which  the  rein- 
forcement was  attached. 
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The  comparative  results  are  given  in  Table  26.1*  The 
dimensions  of  the  slabs  and  reinforcements  and  the  numbering  of 
the  specimens  are  taken  from  [21],  using  the  notation  of  Fig*  25*4. 
Experimental  results  are  taken  from  Table  4 of  [21].  For  two  un- 
cut slabs  the  general  yielding  load  was  found  to  be  380,000  and 
390,000  pounds*  Since  the  lower  value  corresponded  to  a higher 
value  of  the  temperature  than  was  used  in  any  of  the  other  experi- 
ments, we  have  here  used  390,000  pounds  as  the  yield  load.  The 
uniaxial  cutout  factor  \ is  thus  the  yield  load  for  the  slab  in 
question  divided  by  390,000. 

In  each  case  we  have  used  the  "beam11  approach  to  deter- 
mine a lower  bound.  The  lower  bound  for  uniform  edge  load  is 
determined  from  Eqs*  24-.  2 and  24-.  4-,  while  the  lower  bound  for 
uniform  edge  velocity  is  obtained  from  Eq.  25.8.  Finally,  the 
upper  bound  is  obtained  by  considering  sliding  out  of  the  plane 
(Sec,  11)  and  is  given  by 

\ = l-(a+&)+—  • (26,1) 

n 


Lower  Bounds 


Specimen 

Number 

a 

Dimensions  Uniform 
5 Load 

Uniform. 

Velocity 

Experimental 

Value 

Upper 

Bound 

5 

4 

1/72 

1/8 

0,62 

0.74 

0. 83 

0. 86 

6 

4 

1/72 

1/4 

0,62 

0, 74 

0.83 

0.79 

11 

4 

1/4 

1/2 

0.84 

0.87 

0. 92 

1.00 

12 

4 

1/8 

1/2 

0,64 

0.74 

0.85 

0.88 

17 

4 

15/144 

1/2 

0,62 

0. 73 

0»83 

0.85 

18 

i 

3/72 

1A 

0,62 

0.74 

0.87 

0. 88 

Table  26.1 

COMPARISON  OF  RESULTS 
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